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1  2 
Michael  Golomb  and  Joseph  Jerome 

Technical  Sumnary  Report  #  2024 
November  1979 
ABSTRACT 

A  detailed  global  and  local  analysis  of  smooth  solutions  of  the  variational 
problem 

s 

(li)  <5  /  k  (s) ds  =  0  , 

0 

subject  to  position  function  constraints 

(lii)  x(s.)  »  p.,  0  <  sA  <  s,  <••••<  s  <  s  , 

l  i  —  0  1  m  — 

in  2 

is  carried  out.  Here  (p^}Q  c  R  is  prescribed,  x  is  a  vector-valued  func¬ 
tion  with  curvature  tc(s)  at  arc  length  s  and  the  interpolation  nodes  s^ 
are  free.  Problem  (1)  may  be  viewed  as  the  mathematical  formulation  of  the 
draftsman's  technique  of  curve  fitting  by  mechanical  splines. 

Although  most  of  the  basic  equations  satisfied  by  these  nonlinear  spline 
curves  have  been  known  for  a  very  long  time,  calculation  via  elliptic  integral 
functions  has  been  hampered  by  a  lack  of  understanding  concerning  what  precise 
information  must  be  specified  for  the  stable  determination  of  a  smooth,  unique 
interpolant  modelling  the  thin  elastic  beam.  In  this  report,  sharp  character¬ 
izations  are  derived  for  the  extremal  interpolants  as  well  as  structure  theorems 
in  terms  of  inflection  point  modes  which  guarantee  uniqueness  and  well-posedness . 

A  certain  type  of  stability  is  introduced  and  studied  and  Shown  to  be 
related  to  (linearization)  concepts  associated  with  piecewise  cubic  spline 
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ABSTRACT  (continued) 


functions,  which  have  been  studied  for  decades  as  a  simplification  of  the 
nonlinear  spline  curves.  Many  examples  are  introduced  and  studied. 


AMS  (MOS)  Subject  Classifications:  41A0S,  49B15,  49B50,  49F22,  65D10 
Key  Words:  Nonlinear  spline  curves,  elastica,  manifolds  of  extremals, 

perturbation  stability,  mode,  ray  configuration,  rectangular 
conf iguration . 

Work  Unit  Number  6  -  Spline  Functions  and  Approximation  Theory 


SIGNIFICANCE  AND  EXPLANATION 


The  mathematical  formulation  of  curve  fitting  by  mechanical  splines, 
i.e.  thin,  flexible,  elastic  beams  passing  through  freely  rotating  sleeves 
anchored  at  fixed  locations,  is  studied  in  this  report.  These  are  called 
elastica  or  nonlinear  spline  curves. 

As  contrasted  with  the  mathematically  idealized  splines,  which  have 
proven  to  be  of  considerable  utility  and  concerning  which  much  information 
is  available,  the  nonlinear  splines  are  relatively  poorly  understood.  The 
writers  are  attempting  to  understand  and  systematically  construct  these 
curves.  Computer  graphics  obtained  by  other  workers  suggest  remarkable 
efficiency  of  the  elastica  for  curve  fitting.  This  is  perhaps  not  too 
surprising  since  the  nonlinear  spline  represents  an  equilibrium  position 


of  a  thin  beam. 
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EQUILIBRIA  OP  THE  CURVATURE  FUNCTIONAL  AND  MANIFOLDS 
OP  NONLINEAR  INTERPOLATING  SPLINE  CURVES 

1  2 
Michael  Golomb  and  Joseph  Jerome 

fl.  Introduction 

Lot  P  •  {p0.Pi, . . . .p^)  be  an  ordered  set  of  points  in  the  Euclidean  plane  (the 
need  not  be  distinct)  and  let  it  be  required  to  pass  a  smooth  curve  through  these  points  in 
the  prescribed  order.  It  is  an  old  technique  of  draftsmen  to  use  a  mechanical  spline  to 
accomplish  this.  If  the  spline  is  considered  as  a  thin  elastic  beam  of  uniform  cross  section 
with  a  central  fiber  that  is  inextansible,  then  the  strain  energy  of  the  bent  spline  of  length 
s  is  given  by 

s  , 

A  /  tc  (s)ds  +  B 
0 

where  k(s)  is  the  curvature  of  the  fiber  at  arc  length  s  and  A,  B  are  constants.  An 
equilibria  position  of  the  spline  makes  the  energy  functional  stationary,  hence  satisfies 

s  - 

(l.li)  S  I  k  (s)ds  -  0  . 

0 

This  equation  together  with  the  interpolation  conditions 
(l.lii)  x(si)  -  pA,  0  <  sQ  <  Sj  <  •••  <  s#  <_  8 

for  the  position  function  x,  constitute  the  swtheauitical  formulation  of  the  draftsman's 
technique.  The  present  article  deals  with  analytical  (not  graphical  nor  computational)  prob¬ 
lems  arising  from  system  (1.1).  In  elasticity  theory  the  solutions  of  (l.li)  or  of  the  more 
general  equation 
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s  , 

«(/(!♦  x*(e)]ds>  -  0  , 

0 

where  X  is  a  constant,  are  known  as  elastlca.  Their  study  dates  back  to  the  Bernoulli 
brothers,  Euler  and  others  (see  Love  (8,  Ch.  XIX]  for  classical  results).  The  boundary  con¬ 
ditions  in  traditional  elasticity  theory  have  little  in  caaaon  with  the  interpolation  condi¬ 
tions  (l.lii).  To  Materialize  the  latter  ones  in  the  beam  nodel  one  nay  think  or  freely 
rotating  small  sleeves,  anchored  at  the  points  pft,...,p  ,  through  which  the  spline  can  slide 
without  friction.  He  refer  to  the  solutions  x  ■  x  of  the  variational  problem  (l.li,  ii)  as 
extremal  P-interpolants .  In  seme  parts  of  the  present  paper  we  deal  with  extremal  length- 
prescribed  P-interpolants,  in  which  s  -  s^  -  sQ  is  given  in  addition  to  P.  To  materialize 

this  condition  one  replaces  the  sleeves  at  pA,p_  by  pins  which  allow  no  sliding.  In  other 

u  n 

parts  we  consider  extremal  angle-prescribed  P-interpolants,  in  which  the  angles  that  the  spline 

makes  at  p  and  p  with  a  reference  line  are  given.  This  situation  prevails  if  the 

sleeves  at  p.  and  p  are  not  allowed  to  rotate, 
u  m 

Stable  equilibria  positions  in  mechanics  are  sought  as  positions  that  minimize  the 

potential  energy  functional,  in  one  of  the  earliest  papers  discussing  nonlinear  interpolating 

s  _ 

splines  [2]  it  was  pointed  out  that  the  inf  imvxi  of  /  k  (s)ds  is  0  for  any  configuration 

0 

P,  hence  can  be  attained  only  in  the  trivial  case  where  P  is  interpolated  by  a  straight 
segment,  lee  and  Forsythe  [7],  who  make  a  substantial  study  of  the  variational  problem 
(l.li,  ii) ,  call  the  solutions  (when  existence  is  hypothesized)  local  minima.  However,  it  will 
be  proved  in  S6  below  that,  in  the  simple  case  where  P  consists  of  2  points,  there  are 
countably  many  nontrivial  extremal  P-interpolants i  none  of  them  constitutes  a  local  minimum  of 
the  energy.  This  makes  it  evident  that  an  extresMl  P-interpolant  is,  in  general,  not  a  local 
minianim  (for  a  detailed  discussion  of  the  stability  problem  see  [5]). 

The  existence  questions  for  interpolating  elastica  are  much  more  subtle.  For  length- 
constrained  or  length-prescribed  interpolants  one  can  prove  existence  of  extreatals  (actually 
global  minima)  by  the  direct  methods  of  the  calculus  of  variations,  because  one  has  compact¬ 
ness  in  a  suitably  chosen  function  space  (this  was  done  in  [3]  and  (6))  see  also  the 


Appendix  of  this  paper) .  This  Is  not  the  case  for  Interpol ants  with  no  length  restriction, 
and  the  existence  of  such  extremals  interpolating  n  points  in  general  position,  and  whether 
they  are  local  minima  or  not  remains  an  open  question  (same  progress  along  these  lines  has 
been  achieved  by  M.  Golomb  [4],  [5]).  Computational  work  on  extremal  interpolants  is  more 
advanced  (cf.  H.  Malcolm  (10)),  although  decisive  progress  in  this  area  is  al^o  hampered  by 
the  lack  of  general  existence  and  uniqueness  theorems. 

Me  now  give  a  brief  account  of  the  content  of  this  paper.  In  52  we  define  the  function 
classes  in  which  the  extremal  interpolants  are  sought.  He  also  characterise  them  by  Euler 
equations  (for  the  Cartesian  coordinates) ,  boundary  and  regularity  conditions.  In  S3  we  do 
the  same  for  the  "normal  representation”  of  the  extremals ,  by  which  we  mean  the  function 
s  >-*  8(s),  which  is  the  angle  that  the  extremal  makes  at  arc  length  s  with  a  reference  line. 

The  normal  representation  8  of  a  length-prescribed  extremal  P-interpolant  appears  as  the 
solution  of  a  free  multi-point  boundary  value  problem  for  s^ .. . ,Sm-l'  6  with  SB  "  so  pre” 
scribed: 

(i)  8(s)  -  u*  sin  8  (s)  y*  cos  8(s)  -  0,  s^  <  s  <  Sj^  , 

(1.3)  (ii)  8 (s  )  -  8 (a  )  -  0  , 

U  in 

(iii)  (y*+J  -  V^)cos  8(si)  +  (Wx+1  “  V^) sin  BfSj)  »  0,  i  -  l,*”,m-l  . 


For  the  general  extremal  P-interpolant,  (1.3iii)  holds  also  for  i  «  m  and  sA  and  s_  are 

u  in 


free  as  well  with  y 


m+1 


m+1 


0.  The  function  8  and  the  knot  abscissas  si  are  the 


unknowns;  the  multipliers  y^,  y^  are  determined  from  the  interpolation  conditions.  In  54  it 
is  shown  how  certain  families  E..  .  ,  of  extremal  interpolants  with  prescribed  numbers 

see 

(k, ,...,k_)  of  inflection  points  between  knots  ("mode”),  can  be  realised  as  smooth  2m-dlmen- 
x  n 

s Iona 1  manifolds  M,  . .  The  inverse  mapping  from  M..  .  .  into  the  position 

function  space  of  the  extremals  is  continuous  when  the  latter  is  topologised  by  a  suitable 
metric.  Thus,  the  mode  (k^,...,^)  of  an  extremal  suitably  delineates  uniqueness  and  well- 
posedness.  The  union  of  the  E.  consists  of  only  those  extremals  which  have  a 

see 
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genuine  knot  and  no  inflection  point  at  each  interior  interpolation  node.  Points  in  the 

intersection  of  the  boundaries  of  the  manifolds  M  correspond  to  singular  points 

... 

in  position  function  space.  For  these  boundary  elements  some  interior  knot  is  spurious  U  is 
not  discontinuous)  or  is  an  inflection  point.  He  give  two  examples  to  demonstrate  this. 

In  S5  we  study  the  existence  of  elastics  spline  interpolation  in  the  small.  Does  the 
set  of  configurations  P  for  which  extremal  interpolants  exist  have  nonempty  interior  in 
H2"*  ?  More  specifically,  which  P  in  are  interior  points  of  this  set?  We  show,  by 

use  of  the  implicit  function  theorem  that  near  a  given  configuration  P  with  extremal  inter- 
polant  E  there  is  a  local  diffeomorphism  between  configurations  and  extremal  interpolants 
if  (E,P)  satisfies  a  certain  hypothesis  (A) .  It  requires  that  a  homogeneous  linear  differ¬ 
ential  equation  with  variable  coefficients  depending  on  E  and  with  homogeneous  linear  side 
conditions  has  no  nontrivial  solution.  Another  formulation  of  this  condition  is  that  a  com¬ 
putable  function  (involving  many  quadratures)  be  ?  0  at  the  end  point  of  i.  It  is  easily 
verified  that  the  ray  configuration  PQ,  with  the  trivial  extremal  interpolant  EQ,  satis¬ 
fies  (A) ,  so  that  the  existence  of  extremal  P-interpolants  for  all  configurations  P  in 
some  Euclidean  neighborhood  of  any  ray  configuration  is  thereby  demonstrated.  The  differential 
equation  problem  of  hypothesis  (A)  reduces  to  the  natural  cubic  spline  interpolation  problem 
in  the  case  (Pq,Eq) .  This  demonstrates  that  cubic  spline  interpolation  can  be  interpreted 
as  the  result  of  linearization  of  extremal  interpolation  (in  the  sense  of  making  Jk2 ds 
stationary)  near  the  trivial  interpolant  for  the  ray  configuration.  This  proof  makes  precise 
the  old  idea  that  cubic  splines  are  in  some  sense  the  "smoothest”  interpolants.  Of  course, 
it  has  long  been  known  that  cubic  spline  functions  arise  from  minimizing  the  quadratic  func- 

.  2  2  2 

tional  j (D  f)  among  the  interpolating  functions  f .  Since  the  linear  operator  D 
supposedly  approximates  the  nonlinear  curvature  operator,  the  cubic  splines  recommend  them¬ 
selves  as  near  optimally  smooth  interpolants.  The  "hairpin”  configuration  P  with  a  loop 
interpolant  E  is  given  as  an  example  where  hypothesis  (A)  is  not  satisfied.  There  are  con¬ 
figurations  close  to  P  for  which  there  exists  no  extremal  interpolant  near  E  and  there 
are  other  configurations  close  to  P  for  which  there  do  exist  extremal  interpolants 
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near  £.  This  seems  to  be  the  first  known  example  demonstrating  singular  behavior  in  non¬ 
linear  spline  interpolation. 

56  contains  an  exhaustive  study  of  extremal  p-interpolants  for  the  case  where  P  consist 
of  two  points.  It  is  shown  that  there  exist,  besides  the  trivial  extremal,  countably  many  non 
trivial  ones  of  distinct  integral  mode,  that  all  of  them  are  obtained  by  simple  transformation 
from  a  basic  one,  all  have  the  same  length  and  (cf.  [5])  none  makes  the  potential  energy  a 
local  minimum.  Composition  of  these  2-point  extremals  yields  countably  many  extremal  P-inter¬ 
polants  for  various  special  configurations  P.  56  also  exhibits  countably  many  angle- 
prescribed  and  countably  many  length-prescribed  2-point  extremals. 

In  57  some  special  cases  of  closed  extremal  p-interpolants  are  considered.  It  is  shown 
that  the  only  closed  length-prescribed  extremal  without  knots  are  the  repeatedly  traversed 
circle  and  figure  eight  configurations.  Formally,  the  Euler  equation  is  the  limiting  case  cf 
the  Euler  equation  for  an  elastic  circular  ring  under  hydrostatic  pressure  p  as  p  -*•  0 
(cf.  (11  and  (121).  The  ring,  however,  is  not  an  elastics  since  its  deformations  satisfy 
stress-strain  relationships.  We  also  consider  closed  extremals  which  are  not  length-pre¬ 
scribed.  Here  the  extremals  in  56  are  used  to  construct  infinitely  many  closed  extremals 
for  several  special  P-configurations,  for  example  where  P  is  the  set  of  vertices  of  a 
regular  polygon.  In  particular,  if  a  regular  m-gon,  m  >_  2,  is  inscribed  in  the  unit  circle, 
then  a  circumscribed  extremal  exists  with  length 

2m  sin  £  F(|  0^ 

s  -  - : - : -  t 

/2[2E(i  /2;  0a)  -  F<±  /2,  BJ 1 

where  cos  8  »  /cos  it/m.  There  are  similar  formulas  for  the  energy  U  and  the  arc  length 

a  (6);  of  interest  is  the  result  that  s_<8)/8  \  as  m  -*■  •,  so  that  the  circumscribed 

m  m 

extremals  have  the  unit  circle  as  a  limiting  configuration.  These  extremals  are  stable, 
i.e.  they  make  the  potential  energy  a  local  minimum,  as  proved  in  (5] . 


S2. 


ii  nn  j  i-,n^ , .  i.i  5  _-ii  j  .i.ij.  iiii.lm  uyj 


Bsguisrlty  and  Characterisations  of  Open  Extremals 

12  12  2 

For  two  point*  p  -  (p  ,p  )  and  q  -  (q  ,q  )  in  real  Euclidean  space  It  we  employ 

the  inner  product  pq  -  piq1  +  p2q2,  the  distance  |p— q|  “  [  (p-q)  (p-q)  ]  ly^2,  and  the  exterior 

12  2  1  12 
product  [p<q]  ■  p  q  -  p  q  of  such  points.  He  consider  mappings  x  •  (x  ,x  )  of  the  unit 

interval  I  »  [0,1]  to  K2  .  He  denote  by  Hj (I)  the  real  Hilbert  space  of  those  mappings 

x  such  that  the  derivative  x  is  absolutely  continuous  and  x  «  L2 (I) ,  equipped  with  the 

inner  product 

(2.1)  (x,y)  -  /  (xy  +  xy  +  xy)  . 


He  say  x  is  a  regular  element  of  H2  if  |x(t) |  >  0  for  all  t  c  I.  He  observe  that  the 
regular  elements  of  Hj  form  an  open  subset  Hj*9  of  Hj. 

For  x  t  Hj  we  define  the  arc  length  map  i  I  ♦  by 


t 

(2.2)  a  (t)  ■  /  |x|,  t  c  I  . 

0 

If  x  e  Hj"9  then  has  an  inverse  s^1  :  [0,s]  -*■  [0,1],  where  s  -  s^ll) ,  and  in  this 

-1-2 

case  the  function  x  •  :  10, si  -►  B  has  an  absolutely  continuous  derivative,  and  square- 

1  2 

integrable  second  derivative.  He  identify  x  with  the  oriented  curve  c  in  the  x  x  -plane 
which  has  parametric  representation  x  ■  x(t).  Writing  x  ■  x  *  s^1,  we  say  that  x  is 
the  arc  length  parametrisation  of  the  curve  C.  Clearly  x  t  H^CO.s)  and  we  have: 


„  ..  -li»  -1 1  •  -l,’..  -1, 

•  -1  -l.2  -  *•*_  l*»S  I  -  x.s  (xx.s  ) 

X  .  S  /  X  .  s  ,  x  -  - - - - - - — 

x  1  x  1  ,4 


x  .  8_ 


-ll 


Suppose  x  £  and  sx(l)  “  s.  Then  we  define  the  curvature  functional. 


(2.4) 


s  s  •  ..  s  .. 

U(X)  -  /  <2  -  /  1X,X]2  -  /  X  2 

0  0  0 


Note  that  (2.4)  defines  U  as  a  mapping  of  H*®9  into  .  The  equivalent  expression. 


(2.5)  U(x)  »  /  [x,x)2|x|-5  , 

I 

00 

is  independent  of  the  parametrization  of  x  in  the  following  sense.  If  u  is  a  C  -map 
of  I  onto  itself  with  u  >  0  and  *  »  yiu  then 


U(x)  »  /  Cy ,y]  2 ( y  |  ~5  ■=  U(y)  . 
1 


req 

If  x  £  H2  then  U  is  Frechet-differentiable  at  x  and,  for  any  increment  y  £  H2, 

U'(*)[yl  “  /  {2tx,xJ  ( [x,y]  +  (y.xDlxl-5  -  5lx,x]  2xy|x|”7}  . 

I 


(2.6) 


If  the  variable  of  integration  is  chosen  to  be  s  ,  (2.6)  simplifies  to 


(2.7) 


u  u  2-  i  _ 

U'  (x)  ly]  »  /  (2x  y  -  3ieJ|x  y)dsx;  x  ■  x  •  sx>  y  <=  y  •  sx  . 


2 

Let  0  *  p.,p,  be  fixed  points  in  F  ,  not  necessarily  distinct,  but 

oi  m 

p,  ,  /  p, ,  i  «  ],..., m,  and  let  P  denote  the  ordered  set  {p.,p, ,...,p  }.  We  refer  to 
i“i  i  o  l  m 

P  as  a  configuration  in  *2  .  If  x  £  H*89  is  such  that  x(t^)  -  p^^  (i  «  0,1,  ...,m)  for 
some  0£tg<t^<*,,<tl#£l  we  say  the  curve  x  is  an  admissible  P-interpolant.  with 
Knots  pA  (i  «  0,...,m).  The  terminals  x(0),  x(l)  may  or  may  not  be  coincident  with  the 
terminal  knots  p^p^.  The  P-interpolants  defined  here  are  to  be  considered  as  open  even  if 
x(0)  *  x(l).  In  the  physical  interpretation  P^"Pj  for  some  i  ?  j  means  that  the  beam  is 
constrained  to  pass  through  two  sleeves  which  are  fixed  at  the  same  point  p^  but  can  rotate 
independently  of  each  other. 
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Suppose  x  is  a  fixed  admissible  P-interpolant  and  x  «  xo  s^1  is  its  arc  length  para¬ 
metrization,  s^d)  *  *  its  length,  xfs^  -  pi  (i  -  0,1,..., m)  its  knots.  Given  any 
z  e  H2  let  z  m  z  a  ,  be  the  parametrization  of  z  which  uses  the  arc  length  of  x  as 
the  parameter,  and  assume  zia^)  =  0  (i  =  0,1,..., m) .  For  |e|  sufficiently  small,  x  +  z 
is  an  admissible  P-interpolant  and 

(2.8)  U(x+ez)  -  U(x)  »  c  U'(x)[z]  +  o(e)  as  e  •*  0  . 

This  justifies  the  following 

Definition  2.1.  The  admissible  P-interpolant  x,  with  arc  length  parametrization 
x  ”  x  o  s^1,  knots  p^  -  xfs^  (i  -  0,1, length  s  «  s^d)  ,  is  an  extremal  P-inter- 
polant  if 

(2.9)  U*  (x>  [z]  =  0  , 

i.e. , 

,s  u  u  2-  -  -  -1 

J  (2x  z  -  3k  x  z  )  »  0,  z  =  z  o  s  , 

0  x  x 

for  every  z  e  H2  satisfying  zts^  -  0  (i  »  0,l,...,m). 

The  following  proposition  follows  from  (2.8)  by  the  usual  arguments  of  the  calculus  of 
variations.  It  helps  to  explain  the  interest  in  extremal  P-interpolants . 

Proposition  2.1.  Suppose  the  admissible  P-interpolant  x  minimizes  the  curvature  functional 
0  locally,  i.e. 

U(x)  <  U(y) 

for  every  admissible  P-interpolant  y  in  a  neighborhood  of  x  in  H2 ,  Then  x  is  an 
extremal  P-interpolant. 

The  three  major  propositions  of  this  section  follow.  We  use  the  notation  y.  for  the 

J 

restriction  of  a  map  y  to  the  interval  j. 
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Proposition  2.2.  The  admissible  P-interpolant  x  with  arc  length  parametrization  x,  knot:; 
x(s^)  (i  »  and  length  s,  is  extremal  if  and  only  if  the  conditions 

(i)  x  £  C2[0,s],  x(s)  «  0  for  0  <  s  <  s„  and  s  «  s  <  s 

—  —  u  in  —  — 

(2.10) 

1*  2-  2  - 
(ii)  (2  x  +  ac^x)  (s)  «  c^  £  ®  for  s  £  (s^.s^,  i  «  l,...,m 

hold  with  x.-  -  .  £  C°°(s,  ,  ,s.)  (i  « 

(si_1,si)  i-l  i 

Remark  2.1.  Throughout  the  paper  we  use  the  same  symbol  to  denote  regularity  classes  for 
both  scalar  and  vector  functions. 

Proof:  The  implication  (2.10)  -•)  (2.9)  is  routine  and  follows  upon  decomposing  !0,s]  into 
subintervals  determined  by  the  s^,  dot-multiplying  (2.10ii)  by  z,  integrating  by  parts 
and  stunning;  the  continuity  of  xz  ,  the  equations  zfs^)  «  0  (i  *  0, . . .  ,m)  and  the  equations 
of  (2.10)  easily  yield  u' (x) [z]  ■  0. 

Conversely,  if  (2.9)  holds  then,  selecting  z  c  C  10, s)  with  support  in  ' 


i  fixed,  we  have 


i+1  .. 

j_  (2x  +  F)z  *  0 


where  F  ■  3kt  x.  By  elementary  distribution  theory,  (2x  +  F)  -  »  is  in 

x  ‘V  i+l' 

C  (“i'*!-!-].)  and 


D  (2x  +  F)  -  -  .  -  0  . 

(V8i+i’ 


It  follows  that 


(2x  +  3<  *)  .-  -  .  -  c 

x  (•i»8i+1)  i 


and,  recursively,  x  -  -  e  C  (s.,8...).  To  prove  the  continuity  of  x  at  an  interior 

(*i'"i+r  1  11 

knot  select  u  in  c“{0,s)  with  support  in  satisfying  u(s1>  “  0, 


rm 


u' (s  )  •  i  and  put  z  -  (u,0) .  Then,  from  (2.9)  and  integration  by  parts. 


i+1 


0  -  /  ( (2  x  ♦  3k2x) z)  +  2(1,0)  (x(s^  +  0)  -  x^  -  0)) 


i-1 


Since  the  first  term  equals 

ci(z(si)  -  z(si_1))  +  ci+i{z(®i+i)  "  * 


-1  2  - 

which  is  clearly  zero,  we  conclude  that  (x  )  .-  -  is  in  C  (s  ,s  ) .  A  similar  argument 

( ,  8  )  Ora 

-2  u  m 

works  for  x  . 

If  s.  is  either  s_  «  0  or  s  -  s ,  jumps  are  replaced  by  one-sided  limits  and  one 

1  U  IB 

concludes  x(s„  +  0)  -  x(s  -  0)  ■  0.  Assume  now  s  <5.  One  argues  as  above  that 
u  in  m 


(2  x 


2- 

*  3k  x)  - 
X  (8 ,S) 
m 


and  that  x  is  continuous  at  s  .  We  show  that  c  >0.  Indeed,  select  u  e  C°°[0,s] ,  with 

m  m 

u  =  0  for  0  <  s  <  s  ,  satisfying  u(s)  »  1,  u(s)  ■  0,  and  put  z  «  (u,0) .  Then,  from 
—  — ■  m 

(2.9)  and  integration  by  parts  over  (s  ,s] , 

in 

o  -  c  ( z < s )  -  z(5j)  -  c* 
in  in  m 


2 

and  a  similar  result  holds  for  c  .  Thus  c  “0.  By  (2.15i)  of  Proposition  2.4  to  follow, 

in  m 

we  conclude  that  ie2(s)  »  0  for  s  <  s  <  s.  In  particular,  x(s)  -  0  for  s  <  s  <  s, 
x  in  m 

and,  by  continuity,  for  s  <  s  <  s.  A  similar  proof  holds  if  0  <  a. .  This  completes  the 
proof  of  the  proposition. 

We  introduce  the  following  notation  for  the  jump  of  the  third  derivative  of  the  extremal 
x  at  the  knot 

see  erne  erne 

(2.11)  ^x  »  x  lij  ♦  0)  -  x  (5  -  0),  i  -  0,l,...,m  . 

see  see 

If  x  (iQ  -  0)  and/or  x  (sm  +  0)  are  not  defined,  they  are  to  be  replaced  by  0. 
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The  following  proposition  gives  expressi.':*.!*  :>.•  the  oxltiia.it  v.iluo  li(x)  which  do  not 

involve  quadratures. 


Corollary  2.1.  If  x  is  the  extremal  F-interpolant  of  proposition  2.2,  then 

m-1 

(2.121)  U(x)  «  l  c  <p  i  -  p.)  , 

i-0 

and 

m  m. 

(2.1211)  U(x)  =  -2  l  pi.  X 

i-0  1  1 


Proof.  If  we  dot-multiply  ( 2 . lOii)  by  x  and  integrate  over  [0,s] ,  using  integration  by 
parts,  we  obtain 


m-1 

-2U<i)  +  3U(x)  -  l  ci(pi+1  -  pt)  , 


which  is  ( 2 . 121) .  Now  use  (2.10ii)  at  sA  +  0  and  si  -  0  and  subtract  to  obtain 

see 

2A.  x  «  c.  -  c.  .,  which  holds  for  i  -  0,1,..., m  if  we  define  c  -  c  *  0.  By  (2.121) 
x  x  i“X  *1  ro 

we  have 


U(x) 


m 

l  c^p 

i-0 


i+1 


-  P4) 


m 


-  I  Pi<c 

i-0 


i 


ci-l> 


A  see 

-2  l  p  A  5  , 

i-0 


so  (2.12ii)  is  also  proved. 


Remark  2.2.  Since  U(x>  -  0  only  if  x  is  linear,  it  follows  from  (2.1211)  that  an  extremal 

•»  *** 

P-interplant  x  that  is  not  linear  must  have  a  discontinuity  of  x  at  sobm  of  the  knots 

*N  «M 

(or  else,  x  (pQ  +  0)  +  0  or  x  (p^  -  0)  ?  0) . 

From  an  extremal  x,  as  characterised  in  Proposition  2.2,  one  can  obtain  infinitely 
many  other  extremals  by  shifting  the  terminals  x(0)  and  x(l)  along  the  rays  that  are  tan¬ 
gent  to  x  at  pQ  and  p^.  The  value  of  U(x)  is  not  changed  by  these  variations.  Ne 
wish  to  ignore  these  trivial  portions  of  an  extremal  and  will  for  this  reason  adopt  the 
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following  convention.  If  we  speak  of  an  extremal  P-interpolant  x  with  arc  parametrizaticr. 


x.  knots  p^  -  x(s^)  (i  »  and  length  s  then,  unless  stated  otherwise, 

s-  -  0,  s  -  s,  and  the  terminals  are  p.  ■  x (0)  -  0,  p  -  x(i  )  -  x (s)  . 
v  m  0  m  m 

For  some  applications  one  wishes  to  constrain  the  p-interpolants  further  by  prescribing 

the  length  s  of  the  arc  between  the  terminals.  Let  the  class  of  these  P-interpolants  be 

called  length-prescribed  (they  differ  from  the  ‘length-constrained”  interpolants  of  [3]). 

The  next  definition  deals  with  the  extremals  for  U  in  this  class. 

Definition  2.2.  The  admissible  P-interpolant  x,  with  arc  length  parametrization  x  ■  x  °  s^1, 
knots  p^  -  xfs^  (i  «  0,...,m),  length  s  ■  s^d).  is  a  length-prescribed  extremal  P- 
interpolant  if 

<2.13i)  U*  (x)  [z]  +  XS'(x)[z]  -  0 

for  any  z  t  Hjd)  for  which  z®  s^1^)  ■  0  (i  *  0,...,m)  and  X  c  »  determined  so  that 

<2.13ii)  S(x)  i  -  /  ]xf  -  i  . 

I 

For  these  extremals  we  have  a  characterization  similar  to  that  of  Proposition  2.2;  note  that 
S ' (x)  is  given  by 

S’ (x) ly]  -  /  xy,  for  all  y  «  H  (!) t  x  ■  x  •  s*1,  y  «  y  ®  s"1  . 

0  i  *  x 

Proposition  2.3.  The  admissible  P-interpolant  x  with  arc  length  parametrization  x,  knots 
Pj  “  x(s^)  (i  •  0,1,..., m)  and  length  s  -  is  a  length-prescribed  extremal  if  and  only 
if 


(i) 

x  «  C2[0,s],  x(0)  «  0,  x(s)  *  0  , 

(2.14) 

1 

(ii) 

f  7-  -  2 

(2  x  +  3k  x  -  Xx)  -  -  -  c.  e  *  , 

*  X 

1  m  1  f  e  e  •  pOI  f 

m 

m  ess 

(iii) 

Xs  -  U(x)  -  J  c.(p  -  P.  .)  ■  U(x)  + 
i-1  1  1 

2  I  PiA i*  • 

i*0 

I  - 
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hold  with  x (-  -  .  e  C  (s.  ,  ,s.)  (i  -  l»...«m) 

(s^_^,s^)  i-l  i 

Proof.  To  prove  the  implication  (2.14i,ii)  (2.13i)  one  proceeds  as  in  the  first  part  of 

the  proof  of  Proposition  2.2.  If,  next,  (2.14ii)  is  dot-aultiplied  by  x  and  integrated 
over  (O.s)  one  obtains,  using  integration  by  parts, 

n  m  ~ 

-2U(x)  +  3U(x)  -  Is  -  l  c  <p  -  p.  J  “  _2  I  P,A.  x 
i-l  11  i-0  1 

which  give  (2.14iii) ,  which  is  seen  to  be  equivalent  to  (2.13ii)  . 

Conversely,  if  (2.13i)  holds  then  selecting  s  £  c"*[0,sj  with  support  in  tsj/si+i' ' 
i  fixed,  we  have 

-  i-. 

]_  (2x  +  P  -  Xx]s  -  0 


where  F  is  as  in  the  proof  of  Proposition  2.2.  It  follows  that 


(2*  +3V*  -Ci 


The  regularity  properties  of  x  are  proved  as  before.  The  argument  in  the  first  part  of  the 
proof  shows,  that  given  (2.14i,ii),  then  (2.14iii)  and  (2.13ii)  are  equivalent. 


Remark  2.3.  For  any  configuration  P  -  {pQ,p^, . . . ,pB) ,  there  exists  an  extremal  P-inter 
polant  satisfying  (2.14)  with  X  e  »  .  Indeed,  if 


m-1 

L0  "  l  lp 

i-0 


i+1  -  Pi' 


then  the  length-prescribed  extremal 
polants  with  length  equal  to  L  >  LQ 


f 


x,  which  minimises  /  ic  among  all  admissible  P-inter- 


is  guaranteed  to  exist  [6]  and  satisfies  ( 2 . 13i)  (cf. 


Appendix] . 


-  2 

Remark  2.4.  Curves  x  :  K  -*•  K  ,  satisfying  the  equation 

-•  2-  - 
2x  +  3k  x  -  \x  *  c 
x 

are  called  elastica  (cf.  [8]),  more  specifically  inflexional  elastica  if  the  curve  has  inflec- 

2  2 

tion  points  (which  is  the  case  if  and  only  if  X  <_  c  ) .  Curves  for  which  X  ■  0,  the  case 
of  primary  interest  in  this  paper,  will  be  referred  to  as  simple  elastics.  Geometrically,  sim¬ 
ple  elastica  are  characterized  by  the  property  that  the  angular  variation  between  consecutive 

inflection  points  is  exactly  it  (for  all  inflexional  elastica  the  angular  variation  is  ^  n) . 

2 

A  smooth  oriented  curve  in  S  with  continuous  curvature  which  consists  of  finitely  many 
subarcs  of  the  simple  elastica  and  has  (possibly)  discontinuities  of  the  curvature  derivative 
at  the  interpolation  points  only  is  called  an  interpolating  elastica  (cf.  [8]). 

The  next  proposition  deals  with  implications  and  equivalences  of  (2.14ii) .  It  should 
be  observed  that  these  results  apply  to  Equation  (2.1011)  as  well,  since  the  latter  is  the 
special  case  of  (2.1411)  with  X  -  0. 

Proposition  2.4.  Condition  (2.1411)  implies  each  of  the  following  four  conditions  on 
(s^.s^  (i  -  1, . ..  ,m) : 

2  - 

(i)  <c*  -  c^  +  X  , 

(ii)  ix  -  ctl  , 

(2.15) 

(iii)  <x  -  j[x,  cil  +  Yj,  «  *  , 

(iv)  •'x  +  ^  Kx  “  <V2)kx  -  0  . 

Moreover,  (2.151)  also  implies  (2.1411). 

Proof;  If  (2.1411)  is  dot-multiplied  by  x,  one  obtains,  since  |x|  ■  1,  |x|  »  tc*, 

•••  M  _ 

«  •  0,  **  ♦  1*1  •  °! 

a  •••  •  M  .  *  ^  1 

3k*  ♦  2xx  “k*  +  2|x|*  +  2xx  ■  k*  ■  c^x  ♦  X  , 


-14' 


hence  (2.15i)  .  Differentiating  (2.151)  we  obtain 

•  1  - 
K  X  ■  —  C.X 

xx  2  i 

If  x  (s)  p  0  for  some  s,  then  |x  (s) |  -  |x(s) |  >  0  and  x(s)/x  (s)  is  the  unit  vector 

X  XX 

(-x2(s),  x1  ( s) ) .  Therefore,  x  (s)  ■  -j[x(s) ,  c^J  and  (2.1Sii)  holds  in  this  case.  If 

u.  2.  * 

x  (s)  -0  and  s  is  a  limit  of  s  such  that  x  (s  )  p  0,  then  continuity  of  x  and  x 
x  n  x  n  x 

together  with  the  previous  argument  give  the  same  equation.  On  a  fixed  interval  (s^_^,s.) , 

let  1^  ■  (s  :  <x<s)  I1  0}  «n<i  let  (a,  6)  be  any  subinterval  in  the  deconqsosition  of 

(s.  ,,s.)\f. .  We  must  show  that  (2.15ii)  holds  on  (o,8).  Now  x  and  x  are  zero  on  (a,B). 

1*1  1  1  XX 

e  mm  v  •  ^ 

By  the  continuity  of  x^  on  (s^  ^ , s^)  it  follows  that  x^( B)  »  0  »  (x(8),  c^) .  Moreover, 
since  x^  »  0  on  (a,B)  , 

(2.16)  x(s)  -  as  +  b,  a  €  (a, 6],  a,b  e  »2 

Thus,  from  (2.16), 

[x(s),  c^  «  (a,ci),  s  e  [a, 8)  . 

In  particular,  (a,Cl]  -0  and  thus  (2.15U)  holds  on  (a, 8).  Next,  (2.15iii)  follows  from 

(2.15ii)  by  integration.  To  show  that  (2.15iv)  holds,  let  kMs)  p  0i  then 
-  -7  Ii 

x(s)  ■  x^Js) (-x  (s) ,  x  (s))  and  one  obtains  upon  differentiating  (2.15ii), 

Kx(»)  ■  -  -ix^tsjc^js)  "  ‘  ^(s)l  . 

Thus,  (2 . 15iv)  holds  at  every  point  s  for  which  (s)  p  0.  The  case  when  s  is  a  limit 

point  of  s  for  which  x  (s  )  p  0  then  follows  immediately .  The  case  when  s  is  not  such 
n  x  n 

a  limit  point  is  of  course  trivial. 

To  show  that  (2.15i)  implies  (2.14ii),  assume  x(s)  P  0  for  soon  s  «  (5^,i^+j).  Since 
i  II  ,11  2 

x(a)  end  x(s)/|x(s) |  are  orthogonal  unit  vectors  in  R  ,  we  have,  using  (2.151)  and  the 
fact  that  x(s)x(s)  +  |x(s) )  ■  0, 

1, 

The  authors  thank  S.  D.  Fisher  for  a  helpful  suggestion  here. 
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(c1x(s))x(s)  +  (CjXfx) )x(s)/|x(s) 


(k*(s)  -  X)x(s)  + 


*<*;<«»  M 

- — -  x(s)/fx(s) 


_  e  eee  e  e  •*  ••  ” 

-  (3ic*(s)  -  X) x (s)  +  (2x(s)x(s)  )x(s)  +  (2x(s)x(s)  )x(s)/|x(s)  | 

_  •  «M 

-  (3ic*(s)  -  A)x(s)  +  2x(s)  , 


2 


which  is  just  (2.14ii).  If  s  is  a  limit  point  of  sn  for  which  xfs^)  /  0,  then  the  same 
equation  holds  by  continuity.  If  s  is  not  such  a  limit  point,  then  (2.14ii)  reduces  to 

e  • 

c^  «  -Xx(s),  which  is  clearly  obtained  by  dot-multiplication  of  (2.15i)  by  x(s) .  This  con¬ 
cludes  the  proof  of  Proposition  2.4. 

At  the  end  of  this  section  we  mention  still  another  constraining  condition  for  extremal 
interpolants .  It  consists  in  fixing  the  angle  that  the  interpolant  x  makes  with  a  fixed 
line  at  the  terminal  knot  pA  «  x(s.)  and/or  p_  ■  x(s  ).  Thus,  the  condition  is 

(2.17)  ;<S0)  -  e0  and/or  iiij  -  em 


2 

where  eA,  e  are  unit  vectors  in  *  .  We  refer  to  these  extremals  as  angle-constrained . 

0  IQ 

If  x  is  an  admissible  P-interpolant  with  knots  p^  ■  x(a^)  which  satisfies  the  con¬ 
straint  (2.17)  then  any  other  P-interpolant  in  a  sufficiently  small  H^-neighborhood  of  x, 
satisfying  the  same  condition,  is  of  the  form  x  +  es  where  i  t  Hj, 

*.s“1(si)  -  0  (i  -  0,1,  ...,m),  i»s“1(s0)  -  0  and/or  •  •s"1^)  -  0.  It  is  easily  seen 

that  if  x  is  an  extremal  P-interpolant  with  the  added  constraint  x(sQ)  -  eQ  then  the  “free 
••  • 

boundary"  condition  x(sQ)  ■  0  of  (2.10i)  is  replaced  by  x(*0>  “  eQ.  Similarly  if 
•  _ 

x  ( s  )  »  e_  is  a  constraint  then  this  condition  replaces  x(s )  “  0.  There  is  no  other  change 
in  the  conditions  of  Proposition  2.2. 
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S3.  Normal  Representation  of  Extremals 

In  this  section  the  dependent  variable  is  an  angle.  Let  T(»  T1)  denote  the  l-dir.or.- 
sional  torus.  ♦  e  T  is  represented  by  a  real  number  (also  denoted  as)  $,  one  of  the  stt 
4  +  2k*(k  «  0,  ±1,  ±2,...).  A  continuous  function  6  :  (O.s)  -  T  is  represented  by  a  con¬ 
tinuous  function  (also  denoted  as)  8  :  (O.s)  ■*  R  ,  one  of  the  set  6  +  2kn.  The  derivative 
6  is  always  a  unique  function  (O.s)  -*•  It  .  Let  H^(O.s)  denote  the  class  of  absolutely  con¬ 
tinuous  functions  8  i  (0,5)  ■*  R  for  which  8  e  L,  (0,1).  Then  the  function  x  «  x.  : 

2  D 

(O.s)  ■*  R2  ,  defined  by 

.  s  s 

(3.1)  xUa)  -  /  cos  8.  xf(s)  -  /  sin  8 

6  0  0  0 

-  12 
is  in  (0#s)  and  represents  an  oriented  curve  C  in  the  x  x  -plane,  parametrized  with 

a  • 

respect  to  arc  length,  x(0)  «  0,  x^/x^  -  tan  6,  8(s)  is  the  angle  which  the  curve  C  makes 
1  • 

with  the  x  -axis  at  arc  length  s,  and  8s  is  the  curvature  of  C  at  s.  Conversely, 
given  an  oriented  curve  C  with  cartesian  representation  x  e  H2<0,s),  parametrized  with 
respect  to  arc  length,  there  is  a  unique  function  8^  e  (0,s)  such  that 

x1(s)  «  x1(0)  +  /  cos  8,  x2(s)  -  x2(0)  +  J  sin  8^  . 

0  0 

We  say  that  8^  is  the  normal  representation  (n.r.)  of  C. 

Curves  C  that  differ  by  a  translation  have  the  same  normal  representation.  If  8^  is 

the  n.r.  of  C  then  8^  +  const,  is  the  n.r.  of  a  curve  obtained  from  c  by  a  rotation, 

and  -8  is  the  n.r.  of  a  curve  obtained  from  C  by  a  reflection  at  the  x1-axis.  In  a 

x 

geometric  setting  we  may  identify  curves  C  which  differ  only  by  a  congruence,  and  each  con¬ 
gruence  class  is  represented  by  a  single  function  8  with  the  specification  8(0)  -  0, 

8(0)  >  0  (or  8(0)  >  0  if  8(0)  -0). 

In  many  cases  it  will  be  convenient  to  characterise  extremal  P-interpolants  by  their 
normal  representation.  For  this  purpose  we  replace  Propositions  2.2  and  2.3  by  the  following 
propositions  whose  proofs  we  omit. 
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Proposition  3.1.  The  function  6  £  (0,s)  is  the  normal  representation  of  an  extremal  p- 

interpolant  with  knots  p  (i  »  0,1, at  0  «  5  <  s,  <  •••  <  s  and  length  s  - 

i  u  l  n 

if  and  only  if  the  conditions 


(3.3) 


< i)  6  £  C1[0,s) ,  9(0),  9  (s)  -  0  , 

±2  \  •  2  -  -  - 
(ii)  0  (s)  »  c^cos  8(s)  +  c^in  0(s)  for  s  (s^  ^,s^)  , 

2 

cA  £  3%  ,  i  *  1, *  ,m  , 

s.  •• 

i  ,i 

(iii)  /  cos  0(s)ds  »  p.,  /  sin  8(s)ds  «  p^,  i  »  l,...,m  , 

0  1  0  1 


hold  with  9(-^  -  }  £  C"{5i_1,ii)  (i  » 


Proposition  3.2.  The  function  6  £  H^(0,s)  is  the  normal  representation  of  a  length-pre¬ 
scribed  extremal  P-interpolant  with  knots  p.  (i  »  0,1, ...,m>  at  0  »  s  <  s  <  •••  <  i 

*  0  2  in 

and  length  s  *«  s  if  and  only  if  the  conditions 
n 


(i)  9  £  c  (0,s],  0(0)  -  9  (s)  «  0  , 

*21-2-  _ 

(ii)  8  (s)  *  Cjcos  8(s)  +  c*sin  0(s)  +  A  for  s  £  (s^.Sj) 


(3.6) 


cA  £  *  ,  i  «  1, . . . ,m  , 


r  -  1  1  2 

(iii)  j  cos  0 (s) ds  ■  p*,  /  sin  9(s)ds  »  pf,  i  ■  l,...,m  , 

n  *  a  * 


8  i2 


dv)  a;  -  /  §  -  i  c  (p.-p  >  , 

0  i-1  1  1  1  1 


hold  with 


9<:i-i'V  4  c  (i 


1 1  •  •  •  »w)  • 


I  (ft 


•  •  e  e 

Remark  3.1.  The  conditions  0(s.  -  0)  ■  0(s.  +0)  (i  »  9(0)  =  0,  and  8(s  )  «  0 

i  i  m 

result  in  m  +  1  conditions  on  the  vector  constants  c,,...,c  .  Both  in  Proposition  3.1  and 

1  in 

3.2  we  have 


(3.7) 


11  --  22  -  - 

(i)  (ci+1-c.)cos  9(s.)  +  (c^+^-c^)  sin  -  0  , 

i  =  1, . . . ,m  -  1 

1  2  - 
(ii)  c  cos  6(0)  +  Cj  sin  0  (0)  +  X  =  0  , 

1  -  -  2  -  - 
(iii)  c  cos  0  (s  )  +  c  sin  0(s  )  X  *  0  , 

m  in  in  in 


where  X  is  to  be  taken  as  0  in  the  case  of  Proposition  3.1. 

For  use  in  Sections  4  and  5  wi  state  and  prove 

Proposition  3.3.  The  function  0  t  H1<0.s)  is  the  normal  representation  of  an  extremal 

P-interpolant  with  knots  p.  (i  -  0,...,m)  at  0  •  s  <  »•*  <  s  and  length  s  -  s  if 

i  0mm 

and  only  if  the  conditions 

-  1  -  7  - 

(i)  20 (s)  +  c^  sin  9(s)  -  c^  cos  0(s)  •  0  for 


(3.8) 


s  e  e  » 


11  --  22  -  - 

(ii)  (c^j-c^cos  9(si)  +  (ci+1-ci)sin  0(si)  *  0  , 


(iii)  9(0)  -  0.  0(s)  -  0 


i  *  1#  •  •  •  fin  f 


hold  with  6  -  -  .  «  C  (S.  . ,i.)  (i  «  l,...,m)  and  c,,  »  0. 

i*x  i  m+i 

Proof.  The  forward  implication  foliwws  directly  from  Propositions  2.3  and  2.4.  The  converse 
implication  follows  upon  multiplying  <3.81)  by  0(s)  and  integrating?  if  the  integrated  equa¬ 
tion  is  evaluated  at  s  ■  s^  and  (3.8iifiii)  is  used(  the  constant  of  integration  is  seen  to 
be  0.  Thus,  (3.8)  implies 


i 


j 


i 

> 

i 

i 


i 
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..ju.ipjwwa 


m 


m  ********** 


S2(»)  -  c*  CO*  0(s)  -  c*  ain  6(a)  ■  0  for  s  £  (*i_1,*  ) 

6  £  C1 [0,a] ,  i  -  l....,m  , 

and  th«  rssult  now  follows  froo  Proposition  3.1. 
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§4.  Manifold  of  Extremals 

As  noted  in  the  previous  section,  we  may  consider  equivalence  classes  of  curves  differir. 

by  a  congruence,  with  representer  satisfying  x(0)  -  0,  8(0)  »  0,  6(0)  >  0;  or  6(0)  o  if 

8(0)  -  0. 

Definition  4.1.  The  extremal  interpolant  E  is  proper  if: 

(i)  E  has  nonzero  curvature  at  each  internal  interpolation  node 


(4.1) 


V  1 


1  ,  .  .  .  ,  ®“1  a 


(ii)  Each  internal  interpolation  node  is  a  genuine  knot,  i.e.,  there  is  a  dis¬ 


continuity  y  0  in  the  derivative  of  the  curvature  at  s^,  i  =  l,...,m- 


Definition  4.2.  The  m-tuple  (k  ,  ...,k  )  of  nonnegative  integers  is  the  mode  of  E  if  there 
■  X  m  - 


are  k^  inflection  points  strictly  between  the  (j-l)th  and  jth  interpolation  nodes;  here 


an  inflection  point  denotes  a  point  of  zero  curvature  and  we  note  that  0  must  change  sign. 


For  a  fixed  mode  (k^, . . . .k^) ,  E  ■  E(k  ^  k  }  will  denote  the  class  of  proper  (m+1)- 

1  m 


extremal  interpolants  E  in  the  mode  (k, ,...,k  ). 

1  m 


Proposition  4,1.  E  is  a  (finite-dimensional)  metric  space  under  the  metric 


(4.2) 


d<Bl'V  “  0^1  ldTXl(SE1t)  ‘  dt  *2<ak2t)| 


Here  s  and  x.  represent  the  lengths  and  Cartesian  representations  (parametrized  w.r.t. 

Ej  j 

arc  length)  of  E^ ,  j  *  1,2,  respectively,  where  x^(0)  -  0,  j  »  1,2. 

Remark  4.1.  He  omit  the  routine  proof.  He  observe  that  E  is  not  complete.  Indeed,  if 
each  E.  .  (m  fixed)  is  embedded  in  the  space  of  all  extremal  interpolants,  with 

\X- *••• t K  1 
x  m 

metric  described  by  (4.2),  then  the  boundary  3E  of  E  may  contain  an  extremal  interpolant 
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then  the  configuration 


which  is  not  proper.  We  also  observe  that  if  is  close  to  E 2> 

interpolated  by  is  close  to  that  interpolated  by  E 2-  This  would  not  be  true  if  E^,  E 2 

were  not  restricted  to  a  class  E  .  (see  Example  4.2  at  the  end  of  this  section) . 

i  in 

Now  let  8  C1[0,s  ]  be  the  normal  representation  of  some  E  €  E.  If 

E  E 

0=s<s,<***<s  -  s„  are  the  interpolation  nodes  of  E,  put  8  (s.)  =  a.,  i  =  0,...,m 

0  1  me  Eli 

and  =  (a,,...,a  )  e  Tm  ^ .  Setting  c.  =  -2y.  we  have  from  Propositions  3.1  and  3.3 

E  1  m-1  l  l 

the  existence  of  a  unique  multiplier  y_  =  (y,,...,y  )  £  CR2)m  such  that,  for  k  *  l,...,m. 

Elm 

(i)  ^  02(s)  +  y£  cos  ®E(s>  ♦  sin  eE(s)  “  °>  <  s  <  \  ' 

(4.3)  (ii)  0  <s)  -  „£  sin  0E ( s)  +  y2  cos  8£(s)  =  0,  s  <  s  <  sfc  ; 

and, 

(iii)  eE<0)  -  0,  0E(sE)  -  0  . 

By  introducing  the  more  convenient  notation 

Bk  “  (Ak'8k}  '  \  *  °'  6k  £  t1'  *  "  ' 

where 

vi  -  a  sin  ek'  uk  *  \  cos  Bk  - 


we  may  rewrite  (4.3)  in  terms  of  the  multipliers  for  k  «  l,...,m: 


(i) 

2  ®E<S> 

+  Ak  sin(8E(s)  -  6k)  -  0, 

s,  _  <  s  <  s. 
k-1  k 

(4.4) 

and. 

(ii) 

eE(s>  ♦ 

cos(0E(s)  -  B^)  -  0, 

sk-i  <  8  <  sk 

(iii)  eE(o)  -  o,  eE(sE)  -  0  . 

Since  8-0  is  not  a  proper  extremal  we  must  have  9  (0+)  +  0.  We  consider  it  as 

E 

part  of  the  definition  of  E  that 
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for  all  E  f.  E.  Geometrically  speaking,  E  contains  only  extremals  which  turn  counterclock¬ 
wise  near  the  initial  point. 

Proposition  4.2.  8_  t  C^lO.S-l  is  the  n.r.  of  an  extremal  E  c  E  .  with  inter- 

- -  E  E  IK.  i  •  .  •  #K  ) 

1  HI 

polation  nodes  at  0  •  s.  <  s  <  •••  <  s  “  s„  if  and  only  if 

u  i.  m  e 

A.  8  satisfies  (4.4i  -  iv)  for  some  A^  >0,  6^  e  T1. 

B.  sin(6i  -  Bi+1>  +  0 

..  •  k.+...+k. 

C.  6  (0+)  >  0,  sgn  8  <s.)  »  (-1)  i  »  l,...,m-l. 

E  El 

Proof :  Since  0E  f  C1  we  have,  by  (4.4i) 

(4. Si)  sin(ai-8i)  -  Ai+1  sintttj-Bj^)  “0,  i  »  l,...,m-l  , 

where  »  9E "  ®E*si  +  i£  anti  only  if 
(4 . 5ii)  Aa  coslOj-B^  -  A^+1  costc^-B^)  “  0  . 

The  two  equations  (4.5i,ii)  are  equivalent  to  (4.5i)  and  sinfBj-B^^)  "  0.  Thus,  the  above 
condition  B.  is  equivalent  to  the  condition  that  each  interpolation  node  of  E  be  a  genuine 
knot  (see  4.1ii).  There  are  k^  inflection  points  of  E  between  the  (i-l)th  and  i-th 
interpolation  nodes  if  and  only  if  0E  changes  sign  k^  times  between  *nd  s^»  i*e- 

ki 

(4.6)  sgn  <i_1  •  sgn  xA  -  (-1)  ,  <1  -  0^8^  • 

Since  by  (4.4iv)  8  (s)  >  0  for  all  sufficiently  small  s,  (4.6)  is  equivalent  to  the  above 
condition  C.  Condition  C.  also  implies  that  ^  0  for  i  *  l,...,m-l,  thus  condition 

(4.1i)  is  also  satisfied. 


Remark  4.2.  Condition  B.  also  implies 

<4‘7>  Ai  -  Ai+l  '  °'  1  *  1 . m-1  • 

Indeed  if  -  Ai+1  -  0  then  by  (4.5i)  sinta^-e.)  •  sinfc^-B^) ,  hence  Bi  =  6  or 
B^  “  ®i+l  +  *  (mod  2ir) ,  which  contradicts  B. 

He  now  define 

K  = {A  £  <  !  Ai  -  Ai+i  *  °-  1  - 1 . ra-i) 

(4.8)  ?  ■  (8  e  T*  i  sin(8^-6i+1j  f>  0,  i  -  l,...,m-l) 

B  -  BE  -  (Blt...,BB).  «  -  -  <V"“m)  • 

By  Proposition  4.2  each  extremal  E  e  E.  determines  a  unique  point  B  e  JRm  * 

±  m 

In  the  next  three  propositions  we  shall  describe  the  mapping  E  ■+■  B  via  the  composition  of 

E 

two  mappings;  the  homeomorphism 

<4.9i)  J  s  E  -*•  (oe,Be> 

of  E  onto  J(E)  c  ^  1  *  (i®  x  Tm)  i  and  the  projection 
(4.9ii)  M  j  (u,B)  -►  B 

of  J(E)  into  i"  «  T®  ,  which  is  a  local  diffeomorphism.  The  composition  M»j  is  a 
global  homeomorphism. 

Proposition  4.3.  The  mapping  J  is  a  homeomorphism  of  E..  onto  its  image. 

Proof;  The  continuity  of  J  follows  directly  from  (4.2)  and  (4.4);  note  that  B,  ,  ...,B 

l  m 

can  be  expressed  via  (4.4)  in  terms  of  a^,  k^,  A^,  thus  also  in  terms  of 

0_  “  arc  tan  x_  ,/x_  Suppose  now  that  J (E)  -  (a_,B_).  We  show  that  E  e  E  is  uniquely 

S  aif  1  E  E 

and  continuously  determined  by  its  map  (a_,B_) .  By  (4.4) 

E  E 

sin(6E(0)  -  8^)  ■  0,  cos(0E(O)  -  81>  <  0  , 
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hence  o„  ■  6_ (0)  ■  g,  +  it  (mod  2ir) .  The  restriction  of  6„  to  [s.,s,]  is  now  uniquely 
0  E  1  E  0  1 

determined  from 


(4.10) 


1 

(i)  6_(s)  -  I-2A.  sin  (9^,(8)  -  6.  )J 2 

E  1  E  1 

(ii)  eE(0)  -  a0  . 


with  uniquely  determined  from 

(iii)  0  (s^)  “  ®E*S'  “  0  for  values  of  s  in  (O.s^ 

Indeed  if  there  is  an  s'  >  s  for  which  9  (s,)  -  9_(s')  *  a  then  9  (s)  *  0  for  some 

<  s  <  s^,  hence  8e(s)  -  0  for  more  than  values  of  s  in  (0,sp  .  This  clearly 

leads  to  an  inductive  process)  indeed  if  0£  is  defined  on  [O.s^l ,  one  obtains  the 

restriction  of  9  to  Is. ,s  1  from  the  initial  value  problem  defined  by  (4.4ii)  with 
E  1  1+1 

initial  values  e^s^)  and  ^(s^) .  a^+1  is  uniquely  determined  from 

8„(«)  "  0  for  k .  values  of  s  in  (s  ,s  )  . 

E  i+1  i+1  E  i+1  x  i+I 


The  process  is  terminated  at  i  -  m-1  by  replacing  the  condition  0E<Bi+1)  *  ai+1  by 
*  -l 

8_(s  )  ■  0.  Since  the  continuity  of  J  is  an  easy  consequence  of  (4.4)  the  proof  is  com- 
E  111 

plete. 

We  determine  now  the  image  set 


(4.11) 


V' 


J(E(k! . v’ 


The  following  are  necessary  conditions  for  (a,B)  e  St 
(i)  sin(a^  -  8^)  <  0  i  ■  l,...,m-l  . 

k,+...+k,_  , 

(4.12)  (ii)  If  k^  ■  0  for  some  2  <  i  <  m-1  then  (-1)  sin(a£  -  a^)  >  0 

(iii)  sin(a^  -  81)  ■  Ai+1  sin(aA  -  Bi+1> .  i  "  l....,m-l  . 
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•  2 

Conditions  (4.9i  and  iii)  express  that  0  is  positive  and  continuous  at  s  #.«.,s  . .  If 

£  X 

some  kj  »  0  then  there  must  be  no  inflection  point  between  end  s^,  hence 

sin(8E  -  a,_^)  does  not  change  sign,  or 

sin(6E<s)  -  a^)  •  ©E CsE_1>  >  0  for  s)[_1  <  s  <_  st  . 

Using  C.  of  Proposition  4.2,  we  obtain  (4.12ii) . 

We  now  show  that  conditions  (4.12)  characterize  the  image  set  S  completely.  We 
observe  that  (4.12i,ii)  define  an  open  set  in  the  (3m-l) -dimensional  space  T™*'1  *  (R®  *  T™) 
while  equations  (4.9iii)  single  out  a  2m-dimensional  surface  in  the  open  set. 

Proposition  4.4.  The  image  set  J(E.V  is 

(JC-  t  •  •  •  »K  ) 

x  in 

S(k  ^  k  .  -  {(a,B)  £  T^-1  X  (*“  x  ?*)  :  conditions  (4.12i,ii,iii>  hold)  . 

1  m 

Proof;  We  need  to  show  that  if  (a,B)  la  such  that  (4.12)  holds,  than  there  are  numbers 
0  »  sQ  <  s^  <  •••  <  s^  ■  sg  and  a  function  0g  e  C*[0,sEl 


such  that  conditions  A.,  C.  of  Proposition  4.2  are  satisfied  and  HK>reover, 
(4.13)  eE<si)  "  ai»  1  “  1*. 


(condition  B.  follows  from  the  definition  of  i^1) .  s^  and  the  restriction  of  8g  to 
[Sq.Sj^]  are  determined  as  in  the  proof  of  Proposition  4.3.  Next,  the  restriction  of  0E 
[s1>s2)  (s2  as  yet  unknown)  is  determined  from  the  initial  value  problem 


0(s)  +  (-1)  1  l~2Aj  ain(6(s)  -  S2)  1 2  -  0 


S (*x)  *  ^  • 


to 


The  solution  is  the  n.r.  of  a  simple  elastics  with  inflection  points  at  equally  spaced 
abscissas  ok  where  6(ok)  *  8j  or  8j  +  *  (mod  2s).  By  (4.121  and  iii)  we  have 
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sinlc^  -  Sj)  <  0,  sin(a2  -  82>  <  0  , 

and,  therefore,  0(s)  attains  the  values  and  a2  exactly  once  between  any  two  consecu¬ 

tive  o^.  Thus  if  there  are  k2  ^  1  inflection  points  between  s^  and  s2  there  is 
exactly  one  s2  for  which  0{s2>  »  a^.  If  k2  “  0  and,  say  Bfs^  >  0,  then  by  condition 
(4.12ii)  sin(a2  -  a^)  >  0,  which  implies  that  0(s)  attains  the  value  a2  for  s2  >  s^, 
with  no  inflection  point  between  and  s2>  By  the  same  arguments  the  values  of 

s.,,...,s  .  and  the  restriction  of  8_  to  [s, ,s , . . . ,  [s  ,,s  ,]  are  determined,  s  «  s„ 

J  m-i  E  3  4  m-2  m-1  m  E 

and  0  on  [s  ,s  ]  are  similarly  obtained,  except  that  the  condition  0(s  )  «  a  is 
£  m- 1  in  fa  in 

•  i 

replaced  by  6(s  )  =  0.  The  obtained  function  8  is  in  C  [0,s  1  because  of  condition 

m  EE 

(4.12iii),  and  it  satisfies  (4.13)  and  conditions  A.  and  C.  of  Proposition  4.2  by  construction. 
Proposition  4.5.  The  projection  M | ^  is  a  local  diffeomorphism  onto  an  open  subset 

k  j  of  x  Tm  .  Thus  S  is  a  2m-dimensional  smooth  (even  analytic)  manifold 

(S  is  not  connected  if  m  >_  1) .  The  composition  map  M  •  J  is  a  (global)  homeomorphism  of 


E(k. . k)  onto  V . k  )• 

lm  1  m 

Proof:  Choose  any  (a°,B°)  t  S,  a°  «  <a°. 


0.  .  „0 
,a  )  and  B 
m 


a0  A0 

'V  81' 


,  8  ) .  Let  U 

ID 


denote  the  open  subset  of  Tm  1  x  (b™  x  im)  satisfying  conditions  (4.12i,ii) .  Define  a 


mapping  •f  :  U 


_m-l 

B  by 


(4.15) 


^(o.B)  »  At  sin(ai  -  8^  -  Ai+1  sin  (cu  -  8i+1>  .  i  “  1 . m-1 


With  this  notation,  (o,B)  e  U  is  in  S  if  and  only  if  ^(a,B)  -  0  (i  ■  l,...,m-l).  Now 
^(a°,B°)  »  0  and  the  Jacobian  [3^/3c^]  is  nonsingular  at  (o^,B°),  since  it  is  a 
diagonal  matrix  with  diagonal  entries 


(4.16) 


^  (a°,B°)  -  aJ  coXaJ  -  ej)  -  Aj+1  cos<a;  -  B°in) 

-  ^  0  . 


We  conclude  there  is,  for  every  neighborhood  c  «  of  (a°,B°)  a  neighborhood  NQ  of 
B°  in  B™  *  t"1  and  a  ^-mapping  a  of  NQ  such  that  o(B°)  -  a°  and  ^(a(B),B)  -  0 
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This  proves  that  M  is  a  local  diffeomorphism. 


for  all  B  f  NQ. 

By  Proposition  4.3  the  mapping  J  is  a  homeooorphism  of  E  onto  S.  If  t  .■ 

(K_  fee*  fK  ) 

i  m 

composite  map  M  « J  is  not  a  homeooorphism,  there  must  be  (a,B),  (a',B)  in  S  with 

a  j  a 1 ,  say  a.  ft  a!  (mod  2s).  By  (4.12iii)  we  have  A.  sin(a.  -  6.)  -  A.  ,  sin(a.  -  i 

11  1X1  1+1  1  1+1 


0  and  A±  sin(a’  -  6^  -  Ai+1  sin(a|  -  Bi+1> 


0.  These  equations  imply  a| 


a .  +  -  ( nod 

l 


2s) .  But  by  (4.12i),  sinta^  -  8^)  <  0  and  sin(a|  -  B^)  <  0,  which  contradicts  the  pre¬ 
vious  conclusion.  Thus  Proposition  4.S  is  completely  proved. 

Remark  4.3.  It  seems  to  be  difficult  to  give  an  intrinsic  characterization  of  the  set 

k  j .  Examples  show  that  it  does  not  coincide  with  *  Tm.  It  certainly  contains 

1  m 

points  B  -  (A, 6)  for  each  combination  of  the  inequalities  A,  <  A_,  A.  <  A„, . . . ,A  ,  <  A  . 

1  2  2  3  m-l  m 

Thus,  M.  and  S  ,  have  at  least  2m  1  disjoint  components, 

in  in 

Ra«— rk  4.4.  The  following  examples  show  that  the  results  of  this  section  fail  if  in  Definition 
4.1  either  (4.1i)  or  (4.1ii)  is  omitted. 


Example  4.1.  Consider  the  3-point  interpolant  eq  with  n.r. 

1 

.2 


$0<s)  -  [-2  sin  8q(s)  ] *  -0,  0  <  s  <  Sj  ■  / 


du 


0  /2  sin  u 


$0(s)  +  (-sin  90<s)  ] 2  -  0,  sx  <.  s  <.  sx  +  / 


du 


/sin  u 


Here  B£  »  (1,0,— ,0),  oQ  »  *,  «  2s,  o2  ■  *.  The  mode  is  (0,0).  EQ  violates  (4.11) 
since  6Q (s^)  -  0.  For  e  >  0  let  the  extremal  E£  of  the  same  mode  (0,0)  be  given  by 
Be  -  (1,0, l,c) ,  so  that 


9  (a)  -  (-2  sin  9  (s)J  ■  0,  0  <  s  <  s,  , 

C  l  *"  *  l|C 


ie(s)  ♦  I-sin{8( (a)  -  c)]2  •  0,  s1>(  <  i  <  s2( 
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If  0  (s,  )  is  close  to  a,  •  2ir  then  0  (s,  )  ■  2tr  -  $  for  some  5  >  0.  Thus,  for 

e  l,c  1  c  l,e 

8  “  8l,e* 

i  1 

0  (s,  )  -  12  sin  6]2  -  -[sin(S  +  e)J2  , 

t  If  b 


which  is  inpossible.  Thus,  a^(B)  cannot  be  defined  as  a  continuous  function  in  a  full 
neighborhood  of  . 

Example  4.2.  Choose  it  <  a„  <  2ir  and  consider  the  3-point  interpolant  Et  with  n.r. 


du 


0,(s)  -  (-2  sin  0  (s))2  -  0,  0  <  s  <.  s  -  /  _ 

*  *  _1*  o  /2  sin  u 

0,(s)  -  (-2  sin  0*(s)l2  -  0,  s,  <  s  <  s  *  2  /  . 

xe  2  o  /2  sin  u 


Here  -  (1,0, 1,0),  aQ  ■  i,  »  a#,  <*2  ■  *.  The  node  is  (0,1).  E,  violates  (4.1ii) 

since  A.  •  A. .  For  *  <  a  <  2s  let  E  be  defined  by  the  n.r.: 

X  « 


0(«)  -  t-2  sin  9(s)3  -  0 1  0  <  s  <  s.  «  /  .  U 

”  “  1  0  /2  sin  u 

•  A  ^  A 

0(s)  -  (-2  sin  0(sn  -  0,  s  <  s  <  i  ,■  2  /  ,  ■ U  . 

1  2  0  /2  sin  u 


Here  B~  -  (1,0, 1,0),  oQ  ■  w,  ■  a,  <>2  -  it,  and  the  mode  is  (0,1)  as  before.  Since  there 
are  extrseuils  for  all  *  <  a  <  2 it,  £  cannot  be  defined  by  B  and  its  mode. 
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.  Perturbations  of  Configurations 


In  the  last  section  it  was  seen  that  the  extremal  interpolants  with  m  variable  inter¬ 
polation  nodes  (more  precisely,  those  that  belong  to  a  fixed  class  E.  )  form  a 

IK- f • • • #K  J 
1  ID 

2m-dimensional  manifold.  One  expects  that  an  arbitrary  configuration  P  «  {0,p, ,...,p  }  can 

1  B 

be  interpolated  by  an  extremal  interpolant  (possibly  by  one  from  each  class  E  ,  .  ) . 

IK_  f  •  •  . 
i  m 

No  solution  of  any  kind  exists  for  this  existence  problem.  In  this  section  we  investigate 
existence  in  the  small.  Does  the  set  of  configurations  P  for  which  extremal  interpolants 
exist  have  nonempty  interior  in  *2m?  More  specifically,  which  P  in  It2111  are  interior 
points  of  this  set? 

To  attack  this  perturbation  problem  one  is  tempted  to  consider  the  rapping  from  the 

2m-dimensional  set  ,  .  that  coordinatizes  the  elements  of  E„  ,  .  (see  5  4), 

(k  , ...,k  )  (k, ,...,k_ ) 

x  in  in 

2n 

or  from  another  2m-dimensional  set  of  parameters,  to  the  configurations  in  R  which  are 
interpolated.  However,  this  mapping  is  so  complicated  -  it  involves  the  elliptic  integrals 
which  are  the  solutions  of  the  extremal  equations  -  that  little  insight  is  gained  from  its 
consideration.  For  this  reason  we  start  with  the  extremals  themselves,  as  defined  by  their 
differential  equations. 

Let  8  :  (0,s  ]  •»  T1  be  the  normal  representation  of  a  given  extremal  interpolant  E, 

m 

which  interpolates  the  configuration 


so  that 

(5.1) 


P  *  <0.p, . 5  }  , 

1  m 


/  cos  8  «  p7,  /  sin  8  •  p  f ,  i  -  0,1,..., m 
0  0 

0  ■  8-  <  s,  <  •••  <  s_  . 
oi  m 


Since  we  consider  only  extremals  E  with  n.r.  8  near  6,  hence  with  knots  s^ 
we  choose  c  >  0,  e  »  j  min(s^  -  si_1)  and  extend  6  to  the  interval  I0,s],  s  - 
by  setting 


near  s^ 
*m  +  *• 


8 (s)  *  8  (■  ) ,  s  <  s  <  a 

n  b 


30' 


We  introduce  two  spaces  of  mappings  from  the  interval  [0,s] : 

NBV  «  space  of  functions  <  :  [0,s]  -*•  R  of  bounded  variation  V(ic)  and  continuous 
from  the  right  with  k(0)  «  x(s  -  0)  “0  and  norm  V(tc) . 

NBV^  ■»  space  of  functions  9  :  10, s]  -»  T1,  which  are  locally  absolutely  continuous 
and  have  derivatives  6  «  NBV,  with  norm  sup|6|  +  V(8)  . 

Both  NBV  and  NBV^  are  B-spaces.  Clearly  5  as  defined  above  is  in  NBV^  and  6  is  in  NBV. 

If  9  -  NBV^  is  the  normal  representation  of  an  extremal  E  which  interpolates  the  con¬ 
figuration  P  »  {0,p  , . . . ,p  }  at  the  nodes  0  »  8-  <  s.  <  •••  <  s  <1  (more  precisely,  we 
l  m  oi  n 

speak  of  the  linear  extension  of  E  to  length  s)  then  the  following  equations  hold  (see 
Proposition  2.4) : 

“*  1*3 

9  (s)  +  -t  6  (s)  «  0,  for  s  <  s  <  s  ,  i  -  l,...,m 
(i)  2  1 

T (s)  »  0,  for  s  <  s  <  s  . 

(5.2)  (ii)  8(s^  -  0)  -  9(5^  "0,  i  »  l....,m  . 

s.  s. 

1  1  1  2 

(iii)  /  cos  9  -  pi#  /  sin  8  «  p.,  i  »  l,...,m  . 

0  0 

It  is  easy  to  show  that  these  equations  characterise  the  interpolant  E  completely. 

We  rewrite  equations  (5.2)  by  using  the  values 

(5.3)  9 (s^)  »  a^,  8(si  +  0)  «  b^,  i  -  1,...,» 

as  parameters  (but  a.  *  0,  b  •  0  always) . 

u  n 
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(i) 


•  1  t  *3 

0 (•)  ♦  2  /  (s-t)0  (t)dt  ■  •i_1  +  i.  s  <  i  «  1, . .  .,r 

*i-l 


0<s)  -  0, 


S_  <  S  <  S  . 
D  ”  — " 


1  -  1  .3 

(5.4)  (ii)  »i  "  *A  i  *  »»i  i>  -  2  f  (si-t)9J(t)dt,  i  -  1 . m  . 

Vl 

“i  #i 

(iii)  /  cos  0  ■>  p3,  /  sin  0  ■  p2,  i  ■  l,...,m  . 

0  0 

2  in 

Equations  (5.4)  define  implicitly  a  napping  G  from  the  space  P  c  (»)  of  configurations 
P  to  the  space  E  of  extremal  interpolants  E.  To  apply  the  Implicit  Function  Theorem  we 
introduce  a  mapping  G  on  the  produce  space  E  «  P  to  NBV  x  »m  *  K2”1  as  follows.  We  set 

8  “  (0)  s,,...,s  ;  a,,...,a  >  bn,...,b  ),  where 

x  n  x  n  w  m-x 

0  e  NBV^j  s^  t  H,  a^  e  R,  b^  £  B 
D  -  NBV^  x  n(5t-I,  Sj+e)  *  *“  *  / 

and  define  a  mapping  G  -  (g,  r^,  q^)  with  components  g  e  NBV,  ri  e  R ,  e  B2  ,  as 
follow : 

s 


(i) 


g(s)  »  6(b)  +  j  /  (s-t)03(t)dt  -  ai_1  -  b^ts-s^) ,  si_1  <_  s  <  s± 

*i-l 


g(s)  -  0(s) , 


s_  <  s  <  B 
in  — ■  — 


1  a  ^  «3 

(5.5)  (ii)  r^  -  -  a^  -  b^U^-a^)  +  j  /  (Sj-t)©  (t)dt 


'i-l 


(iii)  q,  -  J  eos  0  -  p2,  q2  ■  /  sin  0  -  p2  . 

0  1  0  1 

Clearly,  Equations  (5.4)  are  equivalent  to 

G(«,P)  -  0  . 
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In  particular  we  have  G(§,P)  ■  0,  since  we  assume  that  6  is  the  n.r.  of  the  extremal  ir.t 
polant  £  for  the  configuration  P.  We  need  the  Frechet  differential  <^(©,p)[r]  where 


¥  *  (*(  ^.....t^j  ol,...,ant  Bq.  . . .  .6^) 

is  an  increment  to  0.  The  components  of  C^(S,P) [fl  are  denoted  by  g',  r!,  q!.  One  find 
readily 

13 


g'(s)  -  i(s)  +  ~  j  (s-t) k  (t) ( t) dt  -  ~  ici_i(s~si_i)ti.i 


(i) 


(5.6) 


si-l 


“i-l  "  ®i-l*s_si-l)  +  bi-lti-l’  8i-l  -  s  *  si 


g'  (s)  -  4>(s)  ,  s  <  a  <  a 
in 


(ii)  “  “i  ~  “i.!  "  +  <bi_i  “  2  Ai8ici-l)ti-l  ”  +  f  /  (s^tjK2-^  . 


i-l 


8i  8i 

(iii)  (q3) '  -  -/  iji  sin  0  +  cos  8^  <q2)  *  »  /  *  cos  0  +  tj  sin  6i  . 

0  0 


Here  we  have  used  the  notations 

(5.7)  k  -  0,  0A  -  0(8^  ,  »  0(8^  ,  A^s  -  sjL  -  s^  » 

and  the  relation 

8i 

\  S  e3  -  »(sH  +  0)  -  0(si  -  0)  -  bi  l  -  *(s.  -  0)  , 

8i-l 

which  follows  from  (5.2i)  and  (5.3). 

The  continuity  of  G'  near  (8,P)  is  readily  ascertained  from  (5.6) . 
We  can  now  state  the  main  result  of  this  section. 

Theorem  5.1.  G£(8,P)  is  an  isomorphism  of 

HBV1  *  x  x  »m  onto  mv  »  *'  x  (*2)“ 
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if  and  only  if  (6.P)  satisfies  the  following  hypothesis: 
(A)  The  system  for  the  unknown  *  e.  NBV^ : 


(5.8) 


see  3  -2»  «  • 

P  +  —  <  p  =  0  on  i  = 

(i) 

p  *  0  on  (s  ,s)  , 

m 

s . 

i  1 

(ii)  A.iji  +  A.k  /  p  sin(6-6  )  «  0,  i  *  l,...,m  , 

1  1  0  1 


s . 

i 

(iii)  /  \Ji  cos{8-0.)  ■  0,  i  =  l,...,m  , 

0  1 


has  only  the  trivial  solution  p  -  0.  If  (A)  is  satisfied  then  there  are  neighborhoods 

2  m  • 

V_  in  D  and  M  in  (K  )  of  §  and  P  respectively  and  a  diffeomorphism  n, 

0  P 


n  :  N  +  N 
§  P 

such  that  9  e  Nq  defines  an  extremal  P-interpolant  E  for  P  *  H (0) . 

Remark  5.1.  Explanation  of  the  notations  used  above  and  in  the  following: 

®1  “  Ki  “  5*^),  V»i  «  , 

•  •  • 

AjK  «  K(Sj+0)  -  ic (»i— 0) ,  t^p  -  iji (s^+0)  -  *(ii-0)  , 

•  e  e  e 

in  particular  A  ic  »  K(sn*0)  ,  A  <  »  -ic(s  -0),  etc. 
u  u  m  m 

Proof:  We  first  demonstrate  the  injectivity  of  the  bounded  linear  mapping  G^(S,P)  under 
hypothesis  (A) .  Thus,  assume  that  for  seme  ¥: 

C£(6,P)  If}  -  0  . 

Then  we  have  by  (5.6i): 
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(5.9i) 


(5.10) 


'  3  -2-  „  - 

*  2*  *K  “  0  <si-i,si)'  x  “  !*•••»» 


0  on  (s  .  ,  s)  , 


s . 
i 


Cli.-Q)  +  f  i  (5.-t)^  -  i  j 


i-1 


"  Vl  *  6i-lV  +  bi-lVx  "  0  ' 


and  using  (5.6ii) 
(5. Hi) 

Also  by  (5.6i) 
(5.11ii) 


V(si-0)  ~  ai  +  <(si-0)ti  -  0,  i  -  1, . . .  ,m 


Hi  ( s ^+0 )  -  +  biti  »  0  . 


Since  bi  =  k (s^+0) ,  the  last  two  equations  yield 


(5.12) 


+  tj^ic  *  0,  i  -  . 


The  remaining  equations 

(5.13) 

(5.9iii) 


(qh  '  -  (q?)  *  -  0  (see  5.6iii)) 


ti  «  /  <ji  sin(0-0^)  ,  i-1,. 


s . 

i 

/  i|i  cos<0-0  )  -0,  i  -  l,. 
0  1 


are  equivalent  to 


« fin  • 


•  fin  • 


When  (5.13)  is  substituted  in  (5.12),  one  obtains 

.  .  Si 

(5.9ii)  a  *  +  A  k  /  Hi  sin(0-0  >  -  0,  i-1,...,*  . 

0  1 

By  hypothesis  (A)  the  equations  (5.9i,ii,iii) ,  which  coincide  with  Equations  (5.8),  together 
with  Hi  t  NBV  imply  i|i  .  o  on  [0,sl.  <5.9i)  then  implies  H>  -  0  also  on  (s,s).  Then 

*  *  BI 

by  (5.13),  ti  -  0  (i  -  !,...,*>,  and  by  (5.11),  -  0  (i  -  finally,  by  (5.10), 
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(5.141)  also  requires 


(5.1611)  iMi^O)  -  ♦  tt(«i+0)t1  »  hti^O) 

We  therefore  must  hava 

(5.171)  i(!1tO)  -  i(s^-O)  -  ♦  ^h  ♦  uA,  1  ■  1, . . .  ,m 

Tha  last  aquations  (5.14111)  ara  (aaa  S.6ilx) 

“i  “i 

(5.1711)  -  /  4>  sin  5  ♦  t.  cos  5  ■  vj,  /  i  cos  5  ♦  t.  sin  5  • 

0  1  1  1  o  ill 

Tha  general  solution  qr  i  NBV^  of  (5.141)  is  tha  sum  of  a  particular  solution  and  the 

llnaar  combination  of  3m  functions  with  tha  coafflciants  a,,..., a  ,  . 

x  no  w  ll  n> 

Whan  this  \|i  is  substituted  in  (5.171  and  11)  a  nonhoswgeneous  system  of  3m  equations  for 
tha  unknowns  ,  0^ ,  t^  is  obtained.  Tha  homogeneous  part  of  th'  x  system  corresponds  to  Un¬ 
ease  h  “  0,  u^  «  0,  ■  0,  and  it  has  only  tha  trivial  solution  ai  *  0,  0^  •  0,  tj  •  0, 

as  shown  in  tha  first  part  of  tha  proof.  This  demonstrates  tha  surjectivity  of  c^(A,p) 
and  finishas  tha  proof  of  Theorem  5.1. 

Tha  utility  of  this  thaoram  is  illustrated  by  tha  fact  that  It  readily  implies  the  follow¬ 
ing  important  result. 

Corollary  5.2.  Suppose  P  »  (0,  P^' •••»£„)  i*  tha  ray  configuration  p1  -  (s^O)  (i  -  1 . m) 

with  tha  trivial  interpolant  i.  Than  hypothesis  (JO  is  satisfied,  hence  tha  conclusion  of 
Thaoram  5.1  holds. 

Proof i  0  ■  0  in  this  case.  If  wa  put 

x 

(5.18)  s  -  x,  s  •  x.,  s  -  x,  y (x)  -  /  * 

11  0 

than  y(0)  -  0  and  Equations  (5.8)  become 

y(4)  •  0  on  i  ■  l....»m 

y"  ■  0  on  <»t  ,x> 

(5.19) 

y"(xi+0)  -  y"(*i-0)  -  0,  i  -  1 . m 

ytx^)  «0  i  »  1 . m 
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while  y'  e  NBV^  i.e,  y"  e  nbv,  in  particular  y“(0)  »  0.  These  are  exactly  the  equations 
for  a  natural  cubic  spline  that  interpolates  the  points  (x^O)  (i  »  It  follows 

that  y (x)  =  0 ,  and  the  corollary  is  proved. 

Remark  5.2.  We  briefly  draw  a  connection  between  the  above  corollary  and  natural  cubic  spline 
interpolation.  The  mapping  r  from  the  space  of  configurations  P  to  the  space  of  extremal 
interpolants  E  is  implicitly  defined  by  G(S,P)  »  0.  Perturbation  theory  looks  for  a  pair 
8  -  6  +  y,  P  ■  p  +  Z  (Z  ■  {z, _ _ _ z  })  close  to  the  initial  pair  8,P  for  which 

1  BK 

(5.20)  <^(8,P)m  +  g^,(8,p)[z]  -  0  . 

It  is  readily  seen  that  this  is  system  (5.9i,ii,iii)  except  that  (5.9ii)  is  replaced  by 

9i 

(5.21i)  /  *  cosJS-i^)  -  z^  cos  §i  -  sin  8^  . 


Also,  Equation  (5.13)  is  replaced  by 

Si  _  _  .  _  , 

( 5 . 2 lii )  ti  <•  /  i(i  sin(0-9i)  +  zj  cos  0^  +  z*  sin  9^  . 

Now  suppose  the  initial  configuration  P  and  the  interpolant  I  are  the  ray  configuration 
and  trivial  interpolant,  as  in  Corollary  5.2.  Further  suppose  P  is  the  configuration 

P  -  {(x„,y_),  (x, ,y,),...,(x  ,y _)},  x  »  y  -  0 
u  u  ii  nn  u  u 

with  the  -  z*  and  y^  -  z^  small.  Then  (5.21i,ii)  become  (in  the  perturbation 

approximation) 

xi 

(5.22)  /  1>  -  yA,  4  ^  *  xi  , 

moreover 

S  S  B  S  S 

(5.23)  x(s)  -  /  cos  8  -  /  cos  iji  »  s,  y(s)  -  /  sin  8  <*  /  sin  t|i  »  /  <(i  . 

0  0  0  0  0 


38- 


The  perturbed  extremal  E  is  now  the  graph  of  the  function  x  ■*  y(x) ,  which  satisfies  Equa¬ 
tions  (5.19),  except  that  the  last  equation  is  replaced  by  ytx^)  -  y^  Thus  y  is  the 
natural  cubic  spline  that  interpolates  the  data  (x^,y^>  (i  "  0,l,...,m). 

We  have  shown  that  the  cubic  spline  interpolant  is  the  result  of  linerixation  of  extremal 
interpolation  (in  the  sense  of  making  / *2ds  stationary)  near  the  trivial  interpolant  for  the 
ray  configuration. 

Remark  5.3.  The  differential  equation  that  appears  in  hypothesis  (A) 


<5.24i) 


where  <  »  8  is  defined  by  the  equation 


3  2  • 

+  f  K  *  -  0 


<5.24ii) 


a.  1*3 

6  ♦  j  •  0 


can  be  completely  integrated  by  quadratures  alone.  Indeed  ip  «  1  is  clearly  one  integral! 


(5.25i) 


»  * 


is  another;  for  ^  »  0^  and  0 ^  +  (3/2) 0^5  «  0.  The  third  one  is 


(5.25ii) 


*2(s)  -  sk(s)  . 


In  the  special  case  where  k(s)  £  cs,  three  linearly  independent  integrals  ere 


(5.25iii) 


1,  s,  s"  . 


In  the  remainder  of  this  section  we  discuss  the  replacement  of  hypothesis  (A)  by  a 
simpler  condition,  which  requires  only  that  the  value  of  an  explicitly  given  function  (involv¬ 
ing  many  quadratures)  at  i  be  0.  For  this  purpose  we  introduce  a  condition  on  arcs  of 

simple  elastics.  The  arc  of  E  from  s^_^  to  s^  is  said  to  be  ordinary  if  it  is  either 


straight  or 


«i  !  -  (pj  -  pi_l>I<S1_1,ao°s  ®i  +  *41 

♦  (p|  '  p^.j)  I(si_1)asin  -  x^jcos  5^  *  s^minf*^  -  5j>  +  0  . 
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On*  ****  readily  that  only  exceptionally  is  such  an  arc  not  ordinary.  For  example,  the  art 

* 

from  sQ  ■  0  to  s^,  where  xQ  -  0,  *0  ji  0,  is  not  ordinary  if  and  only  if 

-1  -  -2 

sin  -  Pj^  cos  -  0  , 

i.e.,  the  chord  pQp1  is  tangent  to  E  at  p^. 

Theorem  5.3.  Suppose  the  extremal  E  consists  of  ordinary  arcs  only.  Let  W0)  - l- 

be  the  solution  of  system  (5.8),  exclusive  of  the  condition  on  A  (iK  is  explicitly  con- 

n  u 

structed  below).  Then  hypothesis  (A)  of  Theoren  5.1  is  satisfied  if  and  only  if 

8 

*  m 

(5.27)  i0(5B  -  0)  +  -  0)  /  *0  sin(6  -  §m>  +  0  . 

Proofs  Let  4^,  be  the  integrals  of  Equation  if  *  (3/2) i «  0  (case  <c  /  0)  for  which 


(5.28) 


♦i<«i.l>  -  Vi-l  +  0)  *  °'  V'l-1  +  °>  **  1 


Xi<*i-1)  “  Xi(*i-1  +  0)  “  °*  xi(si-i  +  °>  *  1  • 


These  are  line.  ±  combinations  of  1,  ii^ .  One  finds  easily 


(5.29i) 


♦  i  “  Pit(<i_1)2  -  K1_1<  ♦  Ki_1<S  -  *i_1>  K] 

1  i  .  1  .  3 

xi  ”  pi^"2lci-i'ci-i  +  2Ki-i,c  +  J* *1-1*  (s  "  si-i,ic^ 


where  w*  have  used  the  abbreviation 

(5.29ii)  1/Pi  -  2(xi<>1)2  ♦  i(rw>4  . 


W*  construct  ^  successively  on  the  intervals  ^  “  !«•••»*)•  On  (O.s^  we 


have  since  *Q(0)  -  1,  (0)  -  0t 


(5. 30i) 


*0  "  1  +  V0,*l 


with  (iQ  satisfying  condition  <5.9iii)i 


■40- 


(5.3011) 


/  co«(5-e0>  +  $Q(o)  /  ^  cos<e-eo)  -  o  . 

On*  finds  th*  last  integral  to  b*  p^.  Thus,  !i»Q(0)  is  uniquely  determined  from  (5.30ii), 
and  (5.301)  gives  on  (0,8^).  (If  <  ■  0  on  [0,^1  then  one  finds  4<0 (s)  ■ 

1  -  3(s/s1)2).  Assume  has  been  determined  on  [0,si_1l  (i  <  m) .  Then  il/gts^^)  *nd 

$  (s  -  0)  are  known,  hence  f  (*j_^)  can  be  found  from  condition  (5.8ii). 


(5.31) 


Wl’ 


V®i-r0)  +  Ai-i*  J 


•i-1 


* o  »in(8-6i_1) 


Then  if  k  /  0  on  (s^.s^J  we  have 

(5.321)  *Q  -  *0<5i-l>  *  V*i-l)Xi  +  V'i-l^i  ' 


and  to  satisfy  condition  (5.9iii) < 


■i-1  i 

/  *Q  cos(e-S1)  ♦  IVVj*  ♦  ♦o(Bi-l,Xi,OOB(#"V 

0  "i-1 


(5.3211) 


♦  »o(-*i,lt  *i  co*(®*®i)  "  0 


'i-1 


On*  finds,  using  (5.291)  and 

"i  *1 

£  cos  5  -  pj  -  p[_x,  /_  sin  8  -  PA  -  p2_x 

■i-1  "i-1 


that  th*  last  integral  in  (5.3211)  is  Pi«i  t  0.  Thus  !■  uniquely  determined  from 

(5.3211),  and  (5.321)  gives  *0  on  - 

In  th*  oadtted  case  i  -  0  on  ,  (5.32iii)  are  replaced  by 
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(5.33) 


u>  *0<s>  -  ^(‘i-i’  *  ♦  7 


.  .  .  .  .  i  .  -  -  .  j 

(ii)  /  *0  coije-Bi)  +  VVi1  ( Wi*  +  2  ♦o<Bi-i,(Wi) 


♦  6  V‘i-1* “  0  • 


The  conclusions  remain  the  same  as  before. 


With  9  found  on  (0,s  )  there  remains  condition  (5.9ii)  to  be  satisfied] 
o  m 


(5.34) 


*„<»  -0)  ♦  i<;  -0)  /  sin ( 5-5  )  -  0  . 

u  b  m  0  U  ® 


System  (5.6)  has  a  nontrivial  solution  if  and  only  if  the  constructed  integral  satisfies 

(5.34).  This  proves  the  theorem. 

Example  5.1.  To  illustrate  the  utility  of  the  preceding  theorem  consider  the  configuration 

P  -  {(0,0),  (a,0) ,  (a,b)  } 

with  the  extremal  interpolant  E  whose  n.r.  5  is  defined  byi 


6 (s)  -  0 


0  <  s  <  a 


(5.35) 


0(s)  -  (B/b)  [sin  0(s)]1/,2(  a  <  s  <  s  . 


Here  6  and  s  are  the  definite  integrals. 

v 


B  *  /  sin1^2,  s  ■  a  +  (b/B)  /  sin”1^2  . 

0  0 

That  i  does  indeed  interpolate  the  point  (a,b)  follows  from 

■  *  a  *  i 

/  cos  9  ds  ■  /  (1/9)  cos  t  dt-  (b/B )/  sin“1/2cos  -  0  , 

a  0  0 

s  *  .. 

/  sin  5  ds  -  <b/B>  /  sin  •  sin  '  -  b  . 

a  0 
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Using  the  construction  of  the  preceding  theorem,  one  finds  by  straight  forward  computation : 


K*0(s)  -  1  -  3 (s/a)  2,  0  <  a  <  a  . 

♦  2(s)  -  (b/6)  2(s-a)l(s) .  X2(a)  -  2(b/M2f(s)  , 

*0<s)  -  -  2  -  <12/a)<b/B)25(s),  a<s<i  , 

•  ■  -  .2 
*0(a-0)  +  5(a-0)  /  i(.0  sin <8-5  )  -  jjj  ♦  >  0  . 

0 

By  Theorem  5.2  unrestricted  perturbation  of  the  3-point  rectangular  configuration  P  is 
possible. 

Even  for  this  simple  example  we  were  not  able  to  prove  this  result  in  the  more  direct  way, 
by  expressing  the  parameters  of  the  elastica  spline  in  terms  of  the  coordinates  of  the  inter¬ 
polated  configuration. 

Example  5.2.  Let  p  be  the  "hair  pin  configuration" 


P  -  {(0.0),  (0.S),  (0,0)  >,  6-/(2  sin) 1/2  , 

0 

with  the  extremal  interpolant  E  whose  n.r.  §  is  defined  by: 

9(0)  -  0 

8(s)  -  {2  sin  9<s))l/2,  0  <  s  <  o  i  •  /  (2  sin)_l/2 

0 

-  1-2  sin  9(s)] 1/2,  o  <  s  <  2a  . 

•  • 

Then  i(0)  -  k(o)  -  x(2o)  -  0,  -  x(o+0)  -  k(o-0)  -  2.  One  finds  readilyt 

^(s)  -  j  5(s)  ,  Xjls)  -  k(s)  ,  0  <_  s  <.  o  , 

♦2(s)  -  i(s) ,  Xj(s)  -  x(s),  o  <  s  <  2o  . 

( 1  ,  0  <•<  a  , 

V>  "  ( 

U  ♦  26r(s),  o  <  s  <  2a  . 
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Than 


:  " 

*o(2o-0)  +  k(2o-0)  /  (I.  sin(8-2*) 

0 

6  2  o 

-  -24  -  [/  sin  0  +  J  <1  +  24k) sin  §) 

0  o 

m  -24  -  14-4-24)  -  0  . 

By  Theorem  5.3  hypothesis  (A)  is  not  satisfied,  thus  it  cannot  be  concluded  that  the  hair  pin 
configuration  with  the  extremal  interpolant  B  permits  perturbation.  Indeed,  one  can  show 
directly  that  if  P  is  replaced  by  the  perturbed  configuration  P£  -  {(-e,0),  (0,6),  (e,0)] 
there  exists  no  extremal  interpolant  close  to  E  no  matter  how  small  e  /  o  is.  On  the 
other  hand,  if  P  is  replaced  by  £P  »  {(0,0),  (0,4),  (0,c)},  which  is  also  close  to  P, 
then  there  is  an  extremal  interpolant  ^E  near  E,  E  coincides  with  E  for  the  arc  from 
(0,0)  to  (0,4),  the  remaining  arc  is  the  simple  elastics  joining  (0,6)  and  (0,e).  Thus, 
we  have  an  example  of  singular  behavior  taking  place  in  the  perturbation  from  P  to  £P. 
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16.  Special  Cases  of  Open  Extremals 

He  study  in  seme  detail  in  this  section  extremal  P-interpolants  for  some  special  con¬ 
figurations. 

A.  Two  Point  Extremal  Interpolants 

Let  P  be  the  configuration  P  -  {pg.p^} ,  pQ  ■  0.  For  an  extremal  E  with  normal  repre¬ 
sentation  6,6(0)  -  0,  we  have  by  (3.3): 

(i)  0  £  C^[0,5] ,  6(0)  -  6(0)  -  8(5)  -  0  , 

(6.1)  (ii)  §2(s)  «  c1cos  8(s)  +  c2sin  8(s),  0  <  s  <  s  , 

s  _  i  5  _ 

(iii)  /  cos  8  •  p. ,  /  sin  6  -  p. 

0  0  1 


1  2  -  -  2 
(6.1i)  implies  c  »  0  and  c  sin  6(s)  -  0.  Clearly,  either  c  ■  0,  yielding  the  extremal 

-  -  2  2 

8(s)  3  0,  pQ  ■  0,  Pj^  “  (s,0)  ,  or,  c  )rf  0.  In  this  case  we  write  c  «  -2/i.  Differentiation 
of  (6.1ii)  gives  by  Proposition  2.4  the  differential  equation, 

(6.2)  16(s)  +  cos  8(s)  *  0,  0  <  s  <  5  . 

«»  *  •  — 

When  this  equation  is  integrated  over  (0,5)  and  6(0)  •  9(5)  •  0  is  used,  one  obtains  the 
first  of  equations  (6.1iii)  with  p*  -  0.  From  the  equation, 

(£/2)82  -  -sin  8  , 


it  follows  that  -*  £  5  £  0  if  i  >  0  and  0  <  8  <  »  if  i<0.  Since  the  choice  i  <  0 
amounts  to  a  rotation  through  *  of  the  extremal  corresponding  to  £  >  0  or,  to  a  change  in 
orientation  transforming  6  into  -8,  we  may  assists  i  >  0  and  -»  £  8  £  0.  In  this  case, 
Pj  ”  -d>  ■  -|pQ  -  p^|  in  the  second  equation  of  (6.1iii).  He  may  rewrite  (6.1iii): 

5  5 

(6.1iii‘)  J  cos  8  ■  0,  /  sin  8  »  -d 

0  0 

From  (6. Iii),  1(5)  *  0  is  equivalent  to  sin  8(5)  ■  0.  Since  sin  6(s)  £  0  for  0  £  s  £  s 
end  8  is  continuous,  there  are  only  two  possibilities!  8(5)  ■  0  or  8(5)  -  -i. 
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Altogether  we  have  shown  that  (6.1)  nay  be  replaced  by  the  siapler  system 


(6.3) 


i  -2  •  .  _  _ 

(i)  j  e^ts)  -  -sin  8(s),  0  <  s  <  sj  6(0)  -  0,  6(s)  -  0 

or  -it, 
s 

(ii)  /  sin  8  •  -d  . 

0 


If  s  is  interpreted  as  physical  tine,  16  as  the  displacement  along  a  circle  of  radius  i, 
then  (6.3i)  represents  the  pendulum  equation  with  pendulum  length  1,  unit  mass  and  unit  force 
downward,  starting  from  horizontal  position  with  velocity  0  at  s  «  0  and  reaching  velocity 
0  again  at  s  *  s  when  6  -  0  or  -it.  The  pendulum  swings  from  horizontal  position  6-0 
at  s  -  0  through  one  or  more  half-swings  to  horizontal  position  at  s  -  s.  The  kinetic 
analogue  of  the  elastica  equation  was  discovered  by  G.  Kirchhoff ;  see  [6,  p.  399) . 

One  interpretation  of  the  interpolation  condition  (6.3ii)  is  that  the  time  integral  of  the 
kinetic  energy  jU8)2  divided  by  the  maximum  kinetic  energy,  1,  is  the  prescribed  "minimum 
time"  d.  The  length  of  the  pendulum  is  the  main  unknown  of  the  problem. 

The  solution  of  (6.3i)  for  various  values  of  l  can  be  derived  from  the  solutions  of  the 
same  system  for  1-2.  Indeed,  the  trans forma t ion 

(6.4)  6 (s)  -  6 (/57I  •),  s  -  /t/2  s 

converts  (5.3i)  to  the  "normalized  system": 

:  2  -  .  .  . 

(6.5)  8  (s)  «  -sin  8(s),  0  <_  s  <_  sj  6(0)  -  0,  6(s)  -  0  or  -w  . 

The  solution  of  the  pendulum  equation  (6.5)  with  8(s)  -  -w  is  well  known.  It  is 
explicitly  given  by 

(6.6)  6(s>  -  -  j  -  2  arc  sin(2'1/2sn(2“1/2(s  -  |))),  0<s<s  , 

where  2  ^ 2s  is  the  half-period  of  the  Jacobi  function  sn(u)  -  an(ui  2~1//2)  and  arc  sin 
is  the  branch  of  the  inverse  of  sin  with  range  l-  r,  r] .  Por  a  and  U  -  U(6)  -  /*  62 

*  *  n 


we  find  in  terms  of  the  complete  elliptic  integrals  of  the  first  and  second  kind  [13] : 


(6.7) (i) 


I 


dS 


2/2 ! 


dt 


0  ^sin  6  0  /d-t2)  (i-t2, 


2/2  K(2_1/2)  , 


/2) 


and 


(ii) 


0  -  /  /sin  0  d8  -  2S2  /  —  -  2/2[2E(2-1/2)  -  K(2*1/2)  ]  . 


0  /1-t  /2 


The  analytic  continuation  of  §  (also  denoted  as  0)  to  all  of  »  is  given  by  0(s)  »  -0(-s) 
for  s  <  0  and 


(6.8) 


6<s) 


( 


9 (2s-s)  for  s  <  s  <  2s 


l.9(s-2ks)  for  2ks  <_  a  <  2(k+l)s,  k  *  1,2,... 


It  is  seen  that  0  also  solves  (6.5),  the  value  at  the  new  boundary  point  2s  being  0. 

iu ,  *SJ 

In  general  6  -  (k  *  1,2,...)  solves  (6.5),  with  §(ks)  “  -it  or  0  depending  on  whether 

1^1 KS J 

k  is  odd  or  even.  These  are  all  the  solutions  of  (6.5)  with  free  right  end-point. 

We  now  use  the  solution  0  of  the  normalized  problem  to  express  the  general  solution  of 
the  boundary  value  problem  (6.3i)  for  fixed  s  >  0.  It  is  given  by 


(6.9) 


(i)  0(s)  -  0(k  —  s) ,  0  <  s  <  s 


(ii)  l  -  2(s/ks) 2 


where  k  is  any  positive  integer.  To  satisfy  the  remaing  condition  (6.3ii)  we  snist  n>ve 

s  S  ' 

-d  ■  /  sin  0 (s) ds  ■  /  sin  8 (k  *  s)ds  - 
0  0  • 

-  ks  - 

ttt  /  sin  8(s)ds  *  -  0  1 

Ks  s 


hence 

(6.10) 

independent  of  k. 


s  -  (s/U)d  , 


”47' 


He  write  0^  for  the  solution  corresponding  to  k  (■  1,2,...),  for  the  correspond i:.r- 
parameter  value  in  (6.3i)  and  for  the  value  of  the  curvature  functional  for  6^.  Thus, 


(i) 

I  (•>  ”  8(k  j  s) ,  0  ^  s  £  s  »  (s/U)d 

(6.11) 

(ii) 

ik  «  2(d/kU) 2  -  2 (s/ks) 2,  k  -  1,2,... 

(iii) 

0 -  /  02(s)ds  -  k2U2/d  -  k2(s/i)0  . 

K  Q  K 

All  these  solutions  have  the  same  arc  length  s  »  ds/U  *  2. 2d.  Also  S^ts)  -  6^ (ks)  for 

0  <_  s  <_  s/k  and  0^(s)  ■  0^ (ks-s)  for  s/k  <_  s  <_  2s/k,  etc.  Thus  the  curve  represented  by 

0^  consists  of  k  congruent  arcs,  all  similar  (contracted  by  the  factor  1/k)  to  0^ .  The 

curve  whose  normal  representation  is  6^  is  characterised  as  that  arc  of  the  simple  elastics 

whose  endpoints  are  inflection  points  and  which  has  k  -  1  internal  inflection  points  (it 

belongs  to  class  in  the  notation  of  14)  .  For  the  above  physical  interpretation  the 

result  Mans:  If  the  ratio  of  the  time  integral  of  the  kinetic  energy  to  the  maximum  kinetic 

energy  (in  a  motion  from  the  horizontal  to  the  horizontal  position)  is  to  be  the  fixed  number 

.  2 

d  then  the  pendulum  must  have  one  of  the  lengths  -  2(d/kU)  (k  -  1,2,...)  and  the 
pendulum  makes  k  half-swings  in  total  time  s  -  (s/U)d. 

He  summarize  the  results  in 

Proposition  6.1.  There  jure  countably  many  extremal  P-interpolants  interpolating  a  two-point 
configuration  with  |pQ  -  p^|  *  d  >  0,  one  for  each  k  “  0,1,2,...  .  EQ  is  the  trivial  ray 
Interpol ant:  Efc  is  an  arc  of  the  simple  elastics  with  inflection  points  at  the  terminals  and 
k  -  1  internal  inflection  points.  The  normal  representation  0^  of  E^  is  given  by  (6.11), 
where  0  is  the  elliptic  function  (6.6).  Each  of  the  E^  (k  •  1,2,...)  has  the  same  length 
s  *  (s/U)d,  where  s  and  U  are  given  by  (5.7) .  ■  The  extremal  value  of  the  curvature 
functional  for  is  «  k2^  (k  -  0,1,2,...),  where  U1  «  U2/d. 

Remark  6.1.  It  is  shown  in  [5]  that  none  of  the  extremals  E^,  E^,...  provides  a  local 
minimum. 
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B.  Ray  and  Rectangular  Configurations. 


Corollary  6.2.  Suppose  P  “  {0,p.,...,p  )  is  the  "ray"  configuration  where  p.  3  (s.,0) 

■ —  *  1  n  11 

and  0<s,  <  *■*  <  s  .  Countably  many  nontrivial  extremal  P-interpolants  are  obtained  from 
1  in 

the  2-point  extremals  of  Proposition  6.1  as  follows.  Let  xr.  -  ,  represent  any  of  the  nontrivi 

lU / J 

extremals  for  the  configuration  {0,p. }.  Define  x.-  -  .  as  one  of  the  {p »p,)  inter- 

i-  ^l'S2  *  ^ 

£  .  £  . 

polants  of  Proposition  6.1  or  the  negative  of  it  so  that  x(s^  -  0)  +  0).  Continue  in  thi 

way  to  the  intervals  [s.,s.] , . . . , [s  ,,s  ].  The  obtained  curve  represented  by  x  is  an 

z  j  m-i  m 

extremal  P-interpolant. 


Corollary  6.3.  Suppose  P  =  {p^.p^, . . . .p^}  is  the  "rectangular”  configuration  where  the 

angle  between  p^jf^  and  pipi+i  ^  “  1» ...#m-l)  is  either  0  or  i  n/2.  Let  x  represent  a 

P-interpolant  such  that  the  segment  from  pA  to  Pi+1  (i  =  0, . . . ,m  -  1)  is  any  of  the  extre- 

x 

mals  of  Proposition  6.1  (including  the  trivial  one)  and  so  that  x  is  continuous.  In  this  way 
countably  many  extremal  P-interpolants  are  obtained  for  any  rectangular  configuration  P. 

C.  Angle-Constrained  Two  Point  Extremal  Interpolants 
The  bound «u:y  conditions 

(6.12)  5(0)  •  <Pj-P0)  -  Y0.  5(5)  •  (Pj-P0)  -  Yj.  o  <  <  IPj^-PqI 


are  added  to  the  problem  of  Section  A,  replacing  the  eero  curvature  endpoint  conditions. 

Proposition  6.2.  There  are  countably  many  extremal  { pQ , p^ } - interpolants  constrained  by  con¬ 
ditions  (6.12)  if  yQ  -  y^. 

Proof.  Suppose  6  is  the  normal  representation  of  the  sought  extremal.  It  is  easily  seen 
that  one  can  normalise  8  so  that  8  satisfies  equation  (6.3i),  except  that  the  values  of 
8  at  0  and  5  are  not  0  and  but  given  numbers  8^8^  0  >_  80  >_  -  y,  ■»  6Q  or 

0j  •  80  -  «.  It  follows  that  the  solution  of  the  problem  is  a  symmetric  arc  of  the  curve 
(6.91) ,  with  k  even  if  0Q  -  8lf  k  odd  if  ®0  "  ®1  +  *'  8Cale<5  and  moved  80  that  the 
terminals  art  Pq'P^* 
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D.  Regular  Configurations. 

Suppose  P  «  { Pq > P i > • - • / Pffl }  is  a  "regular”  configuration,  by  which  we  mean  each  segment 

PiPx+l  11  *  0 . m  -  1)  is  of  the  length  d  and  makes  the  same  (exterior)  angle  a, 

1  <-  a  2s,  with  the  following  segment  p^+^p^+2-  We  seek  an  angle-constrained  P-interpolant 

x  for  which 

(6.13)  x(0)  •  (p.-p  )  -  x(s)  •  (p -p  )  -  d  sin  f  . 

xu  m  m  m-i  2 

Corollary  6.4.  For  each  regular  configuration  P  *  {pQ ,p1 , . . . there  are  infinitely  many 
extremal  P-interpolants  constrained  by  the  condition  (6.13),  which  consist  of  m  congruent  seg¬ 
ments  . 

Proof.  Let  z  represent  one  of  the  infinitely  many  extremal  (Pq ,p^ ) -interpolants  EQ  of 
length  s^,  constrained  by 

•  • 

2(0)  •  (Pj^-PqJ  ”  2(5^)  •  (Pj-Pq)  “  d  sin  J 

and  such  that  z(0)  x  (p^-pQ)  =  -  z (s^)  x  (p^-p^) .  Let  x  represent  the  uniquely  defined 

P-interpolant  E  that  extends  E^  by  congruent  pieces.  Then  x  has  continuous  slope  and 
curvature  and  satisfies  (6.13);  hence  E  is  one  of  the  sought  extremals. 

E.  Length-Prescribed  Two  Point  Extremal  Interpolants. 

We  assume  now  the  length  s  of  the  extremal  {p^ , p^}-interpolant  E  is  prescribed. 

Suppose  6  is  the  normal  representation  of  E.  with  the  proper  choice  of  the  coordinate  sys¬ 
tem  we  conclude,  using  Proposition  3.2,  that  0  must  satisfy  the  following  conditions  (assuming 
s  >  d) 

(i)  9(0)  -  1(1)  -  0  , 

(6.14)  (ii)  j  02(s)  «  -sin  0(s)  +  X,  0<.s£s,  i  >  0,  X  e  R  , 

5  s 

(iii)  /  cos  9(s)ds  -  0,  /  sin  0(s)ds  «  -d  . 

0  0 


-SO- 


Sot  t>(0)  '  ^(’  ?  i  *0  '  then 


(0.15)  X  -  ain  9q 

In  particular,  -1  <_  X  <■  X,  We  may  again  interpret  Equations  (6.14)  ae  describing  the  motion 
of  a  t<endulum  swinging  from  some  position  0Q  where  8-0  in  fixed  time  a  to  another 

e 

tvs  it  ion  where  0  •  0  and  so  that  the  time  integral  of  the  kinetic  energy  satisfies  a  certain 

condition. 

By  (6.l4ii),  fl(s)  -  o  if  and  only  if  8(s)  is  8Q  or  -it  -  8Q.  The  condition 

s 

/  cos  0-0  follows  from  (6.1411).  Thus,  (6.14)  may  be  replaced  by  the  simpler  system 

0 

(i)  5 2 ( a )  -  sin  6q  -  sin  8(a),  0  <  s  <_  si  5(0)  “  8^  , 

(b.16)  0 ( s)  ■  0Q  or  -it  -  8Q  , 

8 

(ii)  /  sin  8  •  -d  , 

0 

We  make  the  transformation  (6.4)  again,  to  obtain  the  normalised  aystemi 

1 2  -  *  -  « 

(6.17)  0  (a)  •  sin  0Q  -  sin  0(a),  0  s  n  a,  9(0)  -  Bg,  8(a)  "  “*  ~  80  > 

The  solution  is  qivnn  by 

(6.181)  6(a)  -  (5 (si 8q)  »  -  |  -  2  arc  sinlq  sn(2'l/2(s*;/2) t  q2)  ]  , 

q  -  | ain(it/4  +  0q/2>  |  , 


- 1/2  - 

where  2  a  is  the  half-period  ot  the  Jacobi  function  an(ui  q) t 


-s-0„ 


it+0„ 


(6.1811) 


»  ’  ®<9o>  -  -  /  7= 

u  e  /s 


d9 


de 


in  9.  -  sin  8  /sin  8_  ♦  sin  8 

o  o  o 


2  / 


»/2 


d8 


/sin  8-  +  sin  8 


51 


0  to  n/2.  Clearly  3(8#)  »  0  for  some  8^  between  0  and  w/2  (6t  ”  40*)  .  The  ratio 
d(80)/s(6g)  can  be  seen  to  decrease  monotonically  on  [-it/2,84}  with  values  in  the  entire 
interval  (0,1).  Therefore,  for  any  0  <  d  <  5,  there  is  a  unique  eQ  in  I-n/2,6.)  such 
that  (6.21)  holds.  For  this  unique  8Q,  s  »  s ( eQ)  is  determined  according  to  (6.18ii),  ther. 
8k  and  from  (6.20).  Together  with  X  «  sin  8Q,  these  quantities  satisfy  the  original 

system.  U(80)  is  given  by  (6.18iv)  and  the  value  of  U  for  8^  by 


(6.22) 


-  i"8  -2  2  - 

U.s  -  /  6*(s>ds  -  k  - —  U(8  ) 

K  o  *  S  ° 


Altogether  we  have  proved 

Proposition  6 .3  .  There  are  countably  many  extremal  {p^p^-interpolants  of  pre¬ 

scribed  length  s  >  | P0— P^ I •  Their  normal  representations  are  given  explicitly  by  equations 

(6.18)  and  (6.20),  with  the  angle  of  inclination  8Q  at  pQ  determined  from  (6.21).  For  the 

w  2 

value  of  the  functional  U,  the  relation  1^  *  X  Uj  (k  »  1,2,...)  holds. 

Remark  6.2.  The  curves  of  Proposition  6.2  are  subarcs  of  inflexional  elastics  (cf.  Remark  2.4)  . 
For  illustrations  see  (8,  p.  404,  Figures  48-53). 

Remark  6.3.  In  the  beam  interpretation  the  joint  at  p0  exerts  a  force  R  whose  tangential 
component  is  sin  8^/i  and  whose  normal  component  is  cos  8Q/l .  Thus  R  acts  along  the  line 
joining  pQ  to  p^.  The  magnitude  of  the  force  is  1/i.  For  fixed  s  and  d,  the  force 
R^  in  the  mode  £k  has  magnitude  k2Rj.  For  80  <  tt/2  the  tangential  force  on  the  joints 

is  a  pressure,  for  8Q  >  v/2  it  is  a  pull. 
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37.  Examples  of  Closed  Extremals 

Let  P  ■  {p„»p, , ■ • • «P  )  be  a  configuration  as  in  52.  If  x  £  is  such  that 

u  i  ro  « 

(k)  (V) 

x(t.)  »  p.  (  i  »  0,...,m)  for  some  0  <  t„<  •••  <  t  <1,  and  besides  x  (0+)  “  x  (1-) 
ix  0  m 

for  k  «  0,1.2,  we  say  x  represents  an  admissible  closed  P-lnterpolant  with  knots  p^. 

Suppose  x  represents  an  extremal  closed  P-interpolant,  i.e.  x  makes  the  curvature  functional 
U  stationary  in  the  family  of  admissible  closed  P-interpolants.  Then  the  graph  of  x  is  a 
closed  curve  which  has  continuous  curvature  everywhere,  and  the  curvature  is  continuously 
differentiable  at  all  points  other  than  the  knots. 

In  this  paragraph  we  examine  four  classes  of  closed  extremals: 

A.  Closed  extremals  of  prescribed  length  with  no  knots. 

B.  Closed  synmetric  extremals  with  two  knots. 

C.  Closed  extremals  for  rectangular  configurations. 

D.  Closed  extremals  for  regular  polygons. 

A.  If  x  represents  a  closed  extremal  with  no  knots,  of  prescribed  length  s  >  0,  para¬ 
metrised  with  respect  to  arc  length,  k(s)  its  curvature  at  s,  then  one  finds,  as  in 
Proposition  2.3,  that 

(i)  x  £  C**lO,s] ,  x(k)(0)  -  x(k>(s>,  k  -  0,1,2, 

(7.1)  (ii)  2  x  +  3k2x  -  lx  -  c,  c  e  *2  , 

1*5 

(iii)  X  -  r  /  *  • 

8  0 

Conversely  if  x  satisfies  (7.1)  then  x  represents  a  closed  extremal  with  no  knots,  of  pre¬ 
scribed  length  s,  parametrised  by  arc  length.  For  the  normal  representation  9,  where 

-t  s  _  s 

xA(s)  ■  /  cos  8,  x  (s)  ■  /  sin  9,  (7.1)  gives 

0  0 


(7.2) 


s  s 

(i)  6  t  C*[0,s]  ,  /  cos  §  -  /  sin  §  -  0  , 

0  0 

•  e 

9(1)  -  9(0)  +  2kir  (k  -  0,±1,...)»  9(0)  -  9(5)  , 

-21-2- 

(ii)  0  »  c  cos  9  +  c  sin  0  ♦  X  , 


l  *2 

(iii)  X  -  ~  f  6  . 

3  0 


Proposition  7.1.  For  each  k  »  1,2,...  there  exist  exactly  two  closed  extremals  with  no  knots 

of  prescribed  length  s  >  0.  These  are  the  circle  of  radius  s/(2kn)  transversed  k  times 

and  a  contracted  figure  eight  configuration  traversed  k  times. 

Proof.  We  can  omit  (iii)  in  (7.2)  since  it  follows  from  (i)  and  (ii) .  We  write  A  sin(9  +  a) 

1  -  2  -  - 

with  A  >_  0,  a  e  T  for  c  cos  9  +  c  sin  9.  If  9  represents  an  extremal  then  9  -  u 

represents  the  same  extremal  rotated  by  the  angle  a.  Therefore,  (7.2ii)  may  be  replaced  by 

-2  -  2  2  2 

9  »  A  sin  9  +  Wq,  where  u>Q  >  0.  (uQ  ■  0  means  X  -  0,  which  is  impossible  by  (7.2iii)).  we 

may  also  assume  k  «  0,1, ...,  in  (7.2i).  Thus,  (7.2)  is  replaced  by 

(i)  9  e  c”°[0,s],  /  cos  9  ■*  /  sin  9  «  0  , 

0  0 

(7.3)  9(5)  -  0(0)  +  2kir  (k  -  0,1,2, ...)»  9(0)  -  9(5)  , 

(ii)  92  m  a  sin  0  +  w2,  A  >_  0,  u>Q  >  0  . 


Case  1.  A  «  0.  In  this  case  we  may  assume  0Q  ■  0.  Then  5(s)  -  uQs  and  uQs  -  2k* ,  hence 

uQ  «  2kit/s.  For  k  -  1,2, ...  ,9^(8)  »  2kss/s  satisfies  all  conditions.  9^  represents  a 

-  -  2 

circle  of  radius  s/2k*,  traversed  k  times.  The  value  of  0  for  0^  is  (2k s/s)  . 

2 

Case  2.  A  >  0,  uiQ  >  A.  We  may  assume  A  «  1  since  if  9  is  a  solution  of  (7.3)  for  A  >  0, 

-  -  -  -122  2  -12 
of  length  s  >  0,  then  9  defined  by  9(A  s)  is  a  solution  for  A  •  1,  »  A  w^  >  1,  of 

length  5  -  Al/f2a.  It  9  satisfies  (7.3)  then  9(e)  is  uniquely  defined  by 


9(s) 


tt: 


sln  * 
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(7.4) 


k  in  (7.3i)  aunt  be  positive,  and  uQ  is  uniquely  defined  by 


2kx 


dp 


0  /Uq  +  sin  P 


Por  5  defined  in  this  way  we  have,  after  a  change  of  variable. 


/  sin  6  (s)ds  -  /  a— ■ 

0  -kit 


dy> 


sin  p 


klj*  •*»*...  d,  _  J*  . •**..*. 


0  /a u  +  sin  p 


dip]  <  0 


sin  1 1 


which  contradicts  (7.3i).  Thus  no  solution  exists  for  w*  >  a. 

2 

Case  3.  A  >  0,  wQ  <  A.  We  may  again  assume  A  -  1,  thus  uQ  <  1.  Conveniently  replace  r 
by  0  +  it,  and  write  sin  0,  for  id*,  with  0  <  0,  <  Thus  (7.3ii)  is  replaced  by 

(7.3ii')  0^  m  sin  0*  -  sin  5  . 


In  this  case  -x  -  0*  <_  0(s)  <_  0*,  thus  we  must  have  k  »  0  in  (7.3i).  Since  0(s)  cannot 

•  ^ 

be  monotone,  we  must  have  0(a)  *0  and  sin  0(s)  “  sin  0*  for  some  si  it  is  no  restriction 

_  a 

to  assusie  that  this  happens  for  s  ■  0,  and  9(0)  ■  9t,  9(0)  •  0,  As  s  increases  from  0 
to  some  s*,  §(s)  decreases  from  0*  to  -it  -  9*,  when  0(s*)  “  0.  Since  -it  -  0*  ?  0*  the 
curve  cannot  be  closed  yet,  and  as  s  increases  further,  6(s)  increases  up  to  the  value  0,, 
which  is  attained  for  some  s  -  s**,  and  §(s**)  •  0.  If  we  set  O^t)  ■  0(2s*-t),  we  see 
that  01  satisfies  (7.3UM  and  01(sA)  -  I(s.) ,  hence  0X  -  5,  i.e. 

5(s*  +  t)  ■  5(s.  -  t) 

and,  in  particular,  s,*  -  2s*.  Thus  the  curve  obtained  is  symmetric  w.r.t.  the  point  s  ■  0. 
If  we  put  02(s)  ■  -x  -  9(s*-s)  then  wo  see  that  02  satisfies  (7.3ii')  and 
e2(0)  -  5(0)  -  •*,  hence 

5(s*  -  t)  ■  -x  -  5(t)  , 


i 

i 

i 


4 

5 


M 


the  curve  obtained  is  s yens  trie  with  respect  to  the  line  6 

*** 

The  curve  will  be  closed  iff  the  conditions  /  cos  8 


0  and  8  (st/2)  -  -t,/2. 

8.* 

0,  /  sin  8  ■  0  are  satisfied. 
0 


The  first  equation  follows  directly  f rots  the  symetry  of  the  curve. 

***  ®* 

<7.3i')  0-/  sin  8(s)ds  -  4/  sln  * - d  *  ■ 

0  -s/2  /sin  8#  -  sin 


We  are  left  with 

0  . 


Let  the  last  integral  be  denoted  as  X(8A>.  Clearly  1(0)  »  -  j  0  <  0  and 

I(j)  ■  +<«.  Thus,  there  is  a  value  8k  between  0  and  s/2  for  which  (7.3i‘)  holds,  and  it 

is  easily  seen  that  there  is  only  one  such  value  (approximately  8#  «  40s).  With  this  value 

* 

of  8t  we  have  obtained  a  closed  extremal  E  of  length  sAk.  It  is  an  analytic  curve,  crossing 
itself  at  s  *  ~  sM,  and  consists  of  2  congruent  loops,  each  symmetric  w.r.t.  the  same  axis 
(an  illustration  appears  in  [8,  p.  404]  as  an  example  of  an  inflexional  elastica) .  By  proper 
scaling  the  curve  will  have  the  prescribed  length  s.  The  differential  equation  for  the 
normal  representation  of  the  curve  E*  is 

a2^2  -  sin  8#  -  sin  8^  8^(0)  -  8.,  8^^)  »  . 


Thus,  the  inverse  function  8  ♦  s(8)  is  given  by 


s(8)  -  a  / 


JL£_ 


6  /sin  8.  -  sin  *  2 


where  the  constant  a  is  determined  from 

8. 


T  -  • 


A± 


-*/2 


/sin  8#  -  sir.  V 


•  • 

The  other  extr«mals  Ej.By... 
times  with  scale  factor  ...  , 


in  this  sequence  are  obtained  by  traversing  E^  2,3,... 
thus  their  normal  representations  satisfy 


v*> 


e1(ks),  o<«<s  . 


87" 


A.  In  this  case  the  solution  of  (7.3i)  is  monotonically  increasing  in 


Case  4 .  A„ 


0-  ^ 


s  but  does  not  attain  6Q  +  2tt  for  finite  s. 


Remark  7.1.  The  restriction  of  9^(s)  to  [0,—]  represents  a  length-prescribed  extremal 


(length  «  s/2)  interpolating  the  "loop  configuration"  {pQ.P^}  with  pQ  «  p^  This  is  not  a 

*  s 

closed  extremal,  although  it  is  a  closed  curve.  Each  9^(ks)  (k  «  1,2,...;  0  <  s  <  -)  can  also 
be  considered  as  then.r.  of  such  an  interpolating  extremal.  The  curvature  functional  (potent¬ 
ial  energy)  for  this  extremal  is  seen  to  have  the  value 


2 

< ; 


d 


S  -n/ 2  /sin  8#  -  sin 


i  sin  0,  » 


Clearly  this  function  of  s  has  no  stationary  value.  It  follows  that  there  exists  no  extremal 
interpolant  (unconstrained)  for  the  loop  configuration. 


Remark  7.2.  The  length-prescribed  extremal  interpolating  the  loop  configuration  can  also  be 
obtained  as  the  limiting  case  of  the  extremal  of  Sec.  6.E  as  d  -*•  0.  There  it  was  pointed  out 
that  d(64)  »  0  for  0^  satisfying  (see  6.18iii) 

»/2 


0  -  / 


sin  6 


-9^  /sin  8*  +  sin  0 


d6  . 


Clearly,  this  is  the  above  condition  (7.3i')> 

B.  We  turn  to  the  problem  of  closed  extremals  £  with  two  knots.  We  consider  only 


extremals  that  are  symmetric  with  respect  to  the  line  joining  the  two  given  knots.  We  assume 

•1  »2. 


P0  -  (0,0),  p1  -  (0,-d)  with  d  >  0  are  the  knots  and  that  x  -  (x  ,x  )  represents  the 
extremal  I,  parametrised  by  arc  length.  If  a  is  the  length  of  £  and  x(0)  ■  pQ,  then 
i(k)(5)  -  5<k)(0)  for  k  ■  0,1,2,  and  because  of  the  sysawtry,  x(s/2)  •  p^.  Thus,  we  may 


assume 


ilW 


-x1(s-s),  x2(s)  ■  x2(s-s) ,  0  <_  s  <_  s  , 


•2,s_ 


(7.5)  (7.5) 


x*(0)  ■  x1^)  •  x1(s)  ■  Ot  x2( 0)  “  x2(a)  ■  0,  x2(j)  *  -d 


^(O)  «  xl(s)  -  1,  x2(j)  -  ±1;  x2(0>  -  x2<§)  -  x2(s)  -  0 
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If  f  is  the  normal  representation  of  I  then  by  (7.5)  and  Proposition  3.1, 

O  O  O  O  •  • 

(i)  cos  0(s)  *  cos  0 (s-s) f  sin  0(s)  *  -sin  0(s-s)  , 

0  <  s  <  s  , 

s/2  .  1/2  .  s 

(ii)  j  cos  0  ■  0,  J  sin  0  •  -d,  J  cos  •  »  0  , 

0  0  0 

O 

fS 

/  sin  9  ■  0  , 

0 

o 

(7.6)  (iii)  9(0)  »  0,  I(-|)  «  jn,  I(s)  «  2jir»  j  ■  0  or  -1  . 

•  ♦  0  e 

(iv)  9(0)  »  I(|)  -  Its)  -  0  , 

(v)  02(s)  *  A^cos  9(s)  +  A2sin  9(s),  0  £  s  £  s/2  , 

i  •  j  •  °  # 

■  Ajcos  9(s)  +  AjSin  9(s)  ,  s/2  <_  s  <_  s 

s  »  s/2  in  (v)  gives  A^  »  A*  “  substitution  of  (i)  in  (v)  gives  A2  -  -A2  *  -A.  Thus 
(7.6v)  becomes 

( 7 . 6v ■ )  02(s)  «  -A  sin  Its),  0  <  s  <  s/2  , 

«  o  • 

=  +A  sin  9(s),  s/2  <_  s  <_  s  . 

When  the  conditions  1(0)  -  0,  I(s/2)  -  0  or  -ir,  1(0)  -  I(s/2)  -  0  are  taken  together  with 

(7 ,6v' ) ,  it  is  seen  that  I  [0  j/2)  is  one  of  the  functions  9fc  of  16,  with  s  replaced  by 

«  • 

s/2.  The  symmetry  condition  gives  for  9 j*^2  • j : 

Its)  -  -2 jt  -  I(s-s),  s/2  <  «  <  s  . 

Thus  we  have  found  all  solutions  of  system  (7.6)  and  have  proved 

Proposition  1.2.  There  are  countably  many  closed  extremals  E^E^...  with  n*r. 
with  two  knots  which  are  symmetric  with  respect  to  the  line  joining  the  knots.  B^  is  obtained 
from  the  open  extremal  of  Proposition  6.1  by  reflection  at  the  line  joining  the  knots.  Each 
has  the  same  length  2s,  where  a  is  the  length  of  the  open  B^.  The  value  of  the 
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f 


•  2- 

curvature  functional  U  for  the  extremal  is  2k  where  is  the  value  for  the  oter. 

Er 

C.  In  this  section  we  consider  rectangular  configurations  as  defined  in  Corollary  6.2. 

Proposition  7.3.  Let  P  «  (p.,...,p  <P.}  be  a  rectangular  configuration  as  in  Corollary  6.2. 

0  BO 

There  exist  countably  many  closed  extremals  with  knots  at  p.,...,p  . 

o  m 

Proof.  Since  P  is  closed  there  must  be  an  even  number  of  right  angles  between  consecutive 
segments  pipi+i‘  To  connect  p^  to  P^+1  we  use  either  the  trivial  extremal  or  one 

of  the  2-point  open  extremals  of  Proposition  6.1,  with  the  proviso  that  we  switch  from  one  class 
of  extremals  to  the  other  if  the  angle  at  p^  is  a  right  angle,  otherwise  (if  the  angle  is  0) 
no  switch  is  made.  It  is  easy  to  see  that  infinitely  many  closed  P-interpolants  with  con¬ 
tinuous  curvature  everywhere  can  be  obtained  in  this  way. 

D.  Let  p,,...,p  be  the  vertices  of  a  regular  polygon  ordered  as  they  come  when  the 

x  m 

polygon  is  traversed  counter-clockwise.  Define  pi  for  i  >  m  by  periodicity!  p1  ■  Pi_R|' 

Let  P  ,  "  {p, «P, • • • »P,^_i)  for  k  “  1,2,...  .  P  .  is  a  configuration  of  the  kind  that 
m,JC  X  X+K  l+OK  DfR 

Corollary  6.4  applies  to,  and  if  the  construction  used  there  is  applied  to  ?B  one  obtains 
closed  extremals.  Thus  we  have 

Proposition  7.4.  For  each  regular  configuration  P  .  as  described  above  there  are  infinitely 

m,k 

many  closed  extremal  P-interpolants,  each  composed  of  congruent  segments.  For  each  k,  there 

is  precisely  one  such  extremal  whose  intersection  with  the  polygonal  path  connecting  the  points 

e 

of  P_  .  is  precisely  P_  ..  For  k  -  1,  this  axtramal  *_  circumscribes  the  polygon  counter- 

■pK  UfK  m 

clockwise  and  its  representation  x^  satisfies 

<Pi+x  -  Px)  '  ■  lp1+1  -  pjeos  -,  i  -  l,...,m  . 

If  the  polygon  is  inscribed  in  a  unit  circle,  each  of  the  a  ares  of  may  be  expressed 
in  terms  of  the  inverse  of  its  normal  representation: 
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(7.7i) 


s_(9) 


2  sin  l  F(ii/2,iH 

- ®— 2 -  4  -  r  <  e  <  Z  +  r  . 


2E(i/2',0a)  -  F(y/7,B_)  2  *“  2  m 


where  cos  <i  «  /sin  6  and  cos  B_  •  /cos  In  this  case,  the  length  s  and  the  energy 

m  in 

by. 

,  2m  sin  ~  F(4/T,8J 
9  _ ra  2  m 


*  * 

U  of  E  are  given  explicitly  by 


(7.7ii) 

and 

(7.7iii) 

Finally, 

(7.8) 


B  28(4/2,8  )  -  F(4r,BJ 

t  rn  &  n 


,  *  4n 


:12E(|/T, B  )  -  F(i/T,BJ]2  . 


m  .  * 1  2  b  2  m 

sin  — 


s*<0) 

B 


-*■  1  as  b  -*■  »  , 


so  that  the  extremals  have  the  unit  circle  as  limit. 

Proof.  He  sketch  the  verification  of  (7.7).  Starting  from  the  differential  equation 

2»2 

y^Q  “  sin  0  , 

we  obtain 

e 

*  t  dl  IT  IT 

a  (0)  *  Ym  /  — -T-— ,  aM<  0  <  ir-a_,  a_  ■  —  -  —  . 


a  /sin  * 


ra  —  —  nr  m  2 


Since  the  distance  between  adjacent  points  of  P  ,  is  2  cos  a  ,  the  constant  y  is  deter- 

nip  i  n  i« 

mined  from 

«/2 


Y_  -  cos  a  A  /sin  <f  if 

■  01 


,inS 


-  F^.BJ] 
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This  gives  (7.7i)  .  (7,7ii)  is  immediate  and  (7.7iii)  follows  from, 

* 

*  ,s«/2m  .0  »/2 

0  *  2m  /  9  ds  ■  —  /  /sin  ^  . 

m  '  y  ' 

0  'm 


The  calculation  (7.8)  is  routine. 

* 

Remark  7.3.  It  is  shown  in  (51  that  the  extremals  E  are  stable,  i.e.  they  provide  a  local 


minimum  for  the  curvature  functional 


Appendix 

Let  P  =  {p0,...,pm},  be  given  and  let  LQ  =  £  |p.+1~p.|.  We  assume  that  P  is  not 

i-1  1 

collinear . 

Theorem.  For  every  L  >  LQ  there  exists  a  length-prescribed  extremal  P-interpolant  of  length 
L  satisfying  Definition  2.2. 

Proof .  The  existence  of  a  function  x,  parametrized  with  respect  to  arc  length,  for  which 

U(x)  -  min{U(x)  :  x  is  an  admissible  P-interpolant 
of  length  l) 

is  demonstrated  in  (6] .  (The  modification  required  for  the  ordering  of  the  points  in  P  is 
trivial.)  If  x  denotes  the  closed  subspace  of  H2 [0,L]  consisting  of  those  functions 
vanishing  at  the  knots, 

0  <  s„  <  s.  <  •••  <  s  <  L  , 

—  0  1  m  — 

of  x,  i.e.,  at  those  points  s^^  for  which  x(si>  =  p^,  i  -  0,...,m,  define  f  to  be  the 
mapping  of  X  into  R  obtained  by  f(y)  ■  U(x  +  y)  and  let  H  be  the  function  such  that 
H(y)  =  S(x  +  y)  -  L,  where  S  is  the  usual  length  functional  (cf.  < 2 . 13ii) ) .  Clearly, 

f (0)  «  min{f(y)  :  H(y)  *  0}  , 

and  H*(0)  is  surjective,  since  P  is  assumed  non-collinear.  The  result  follows  from  the 
Lagrange  multiplier  rule  (see,  e.g.,  [9,  Theorem  1,  p.  2431). 
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is  carried  out.  Here  {p^}g  c  *  is  prescribed,  x  is  a  vector-valued  function 

with  curvature  Ms)  at  arc  length  s  and  the  interpolation  nodes  s^  are  free. 
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have  been  known  for  a  very  long  time,  calculation  via  elliptic  integral  functions  has 
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mal  interpolants  as  well  as  structure  theorems  in  terms  of  inflection  point  modes 
which  guarantee  uniqueness  and  well-posedness. 

A  certain  type  of  stability  is  introduced  and  studied  and  shown  to  be  related  to 
(linearization)  concepts  associated  with  piecewise  cubic  spline  functions,  which 
have  been  studied  for  decades  as  a  simplication  of  the  nonlinear  spline  curves.  Many 
examples  are  introduced  and  studied. 


